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1.0 SUMMARY 


As a part of the NASA Langley Research Center program to identi- 
fy the important parameters governing sound transmission into 
airplane interiors, and to determine noise control methods, an 
aircraft interior noise prediction model is being developed by 
Bolt Beranek and Newman Inc. (BBN) of Los Angeles. Ihe work 
includes analytical modeling and integration of information and 
technologies needed to understand sound transmission through a 
fuselage wall into an aircraft cabin. A three phase program has 
been defined for accomplishing the goal. 

In the first phase (now concluded), the basic analytical 
modeling of the transmission problem (interaction of the 
structure with the exterior and interior acoustic fields) was 
undertaken and preliminary validation studies were conducted 
using an unpressurized, unstiffened cylinder as a test article. 
Results of that work were presented in Reference [1]. 

The second phase of work, reported in this document, includes 
the development of the general aircraft interior noise model and 
the laying out of the basic master computer program. Validation 
studies are considered using more advanced test articles (one 
being a stiffened cylinder with a floor partition and interior 
trim, i.e., insulation and lining). 

The third phase, now underway, involves completion of the 
analytical models (including propeller excitation) and software 
development with application to an actual (or simulated) 
aircraft fuselage, along with validation tests, comparisons, 
refinements, and documentation of the finalized model and 
software. 

As stated, this report presents the results of the Phase II 
studies. Hie theoretical developments of Phase I that describe 
the interaction of the structure and the interior acoustic field 
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are generalized to include the effects of sidewall insulation 
and lining (trim). The new analysis leads to a transmission 
coefficient that multiplies the previously derived expression 
for the power inflow to the cavity for the case where the trim 
is absent and an additive term giving an increase in the damp- 
ing of the sidewall structure by the trim. Also a more precise 
calculation of the power flow from a structural mode closely 
coupled (in frequency) to an acoustic mode is developed to 
account for the influence of the radiation damping of the struc- 
tural mode by the highly receptive acoustic mode. 


A major addition is the generalization of the geometry of the 
acoustic space to include a floor (partition) . The complex 
cross-section's modal properties are computed using a finite 
difference approach. Appropriate normalizations and use of the 
data for calculating the acoustic/structural coupling terms and 
the cavity loss factors (using predicted wall admittances) are 
also detailed. 

Comparisons of noise reduction predictions with measurements are 
presented for three test articles: 

1) a ring-stringer stiffened cylinder without floor or trim; 
wall thickness of 0. 000508m (0.020 inches). 

2) a 0.00l6m ( 0.064 in.) thick unstiffened cylinder (the Phase 
I test article) modified with a floor partition and lined 
with a 0.0127m (0.50 in.) thick layer of PF-105 fiberglass 
that is covered with a 0. 0000508m (0.002 in.) vinyl 

film. 

3) a 0. 000508m (0.020 in.) ring-stringer stiffened cylinder 
(same as (1) above) with floor partition lined with a 
simulated trim consisting of a 0.0127m (0.50 in.) thick 
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layer of PF-105 fiberglass covered with a 0.00119ni 
(0.0*17 In.) thick layer of lead vinyl weighing 2. *1*1 kg/m2 
(0.50 lb/ft2), with a 0.0127m (0.50 in.) layer of PF-105 
fiberglass on the inside of each end cap and exposed 
directly to the interior. 

r 

To our knowledge no attempt has ever previously been made to 
compute noise reductions for configurations as complex as cases 
(2) and (3) above. 

2.0 INTRODUCTION 

The present study has the specific goal of developing an analyti- 
cal model that can be used to predict the sound levels in a ring- 
stringer stiffened cylinder that has a partition simulating a 
cabin floor and insulation and lining on the inside of the 
cylinder wall simulating a basic cabin sidewall trim. 

Theoretical developments for harmonic (tonal) excitation and for 
excitation by a reverberant acoustic field are given. 

Predictions of the noise reduction for three different test 
articles are compared to measurements for purposes of model 
validation. No calculations are presented for tonal excitation 
in this report (some have already been given in [1]). Tone 
prediction capability is to be brought to a practical level in 
Phase III, when the propeller induced exterior pressure field 
description is to be incorporated into the model and the modal 
forcing functions derived. 



2.1 Analytical Model 


In Reference 1, a basic discussion of the power flow technique 
adopted for this project is given. To a large extent, the 
concepts used for the predictions of concern here have been 
developed in References [2] and [3]. However, in the Phase I 
report, results from [2] and [3] were specialized to include the 
case of harmonic (tonal) excitation. In this report, the 
concepts previously developed in [1], [2] and [3] are extended, 
mainly by including the trim dynamics. Also, although including 
such information does not represent an advance in the concepts of 
Phase I, much more complex structural and acoustic properties are 
considered. Most notable is the use of finite-difference modal 
data for the cavity. Also the structural (modal) properties of 
the orthotropic cylinder (one stiffened by rings and stringers) 
are utilized. 

2.2 Report Organization 

A brief overview of the organization of this report is included 
here. Basically this report presents results of (1) analytical 
derivations, (2) experimental tests, and (3) validation studies. 

In the analytical development, there is to begin with, one funda- 
mental goal, that being to- incorporate the trim dynamics. This 
has to be done for both the low and high frequency models. After 
that has been accomplished, consideration is given to improving 
the precision of calculations of power flow for certain coupled 
acoustic and structural modes by including the effect of radia- 
tion damping of a highly resonant structural mode when closely 
coupled to a highly resonant acoustic mode. Generalization of 
the tonal transmission calculation is then considered; however 
the question of exterior field for the propeller excitation 
remains for the Phase III study. After that the question of 
cavity modes for the cabin configuration (i.e., that formed by 
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the cylinder wall and floor partition) is addressed. The 
normalization and use of the finite difference data is then 
considered. Special problems in the volume stiffness controlled 
regions of the test articles are then discussed (this extreme low 
frequency regime will lie below the regions of concern when full 
scale fuselage structures are considered) . Finally a brief 
overview of the basic master computer program is given. 

The parts of the report concerned with the experimental work 
describe the tests, hardware, and data acquisition and reduction 
techniques. 

Comparisons of predictions and measurements are presented for the 
noise reductions of the three test articles. The statistical 
analysis of the noise reduction data and comparison technique and 
philosophy are consistent with and identical to that of Ref [1]. 

2.3 Program Management 

The work was accomplished in joint effort by BBN/Los Angeles and 
NASA Langley Research Center. The experimental work was done at 
NASA Langley by C. M. Willis and W. H. Mayes. Mr. Mayes acted as 
LaRC technical representative of the contracting officer (TRCO). 
L. D. Pope served as BBN program manager. Most of Section 4 of 
this report was provided by C. M. Willis. A. G. Piersol of BBN 
assisted with the statistical evaluation of the comparison data 
and M. D. Sneddon of BBN mapped the acoustic modal pattern in the 
cylinder with the various floor partitions. 
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3.0 ANALYTICAL MODEL 


The purpose of the first part of this section is to incorporate 
the effects of sidewall trim (insulation and lining) in the 
transmission analysis. In the present case, the interaction of 
the exterior pressure field, fuselage structure, insulation and 
lining, and interior field is of concern. To integrate the trim 
dynamics Into the analysis, the basic expressions that describe 
the response, transmission, and absorption characteristics of the 
various components of the sidewall system must be considered. 


3.1 Trim Dynamics 


The trim is assumed to be representable by a transfer matrix: 



or 



°U ^11 f I 1 ? 

a v An w,J 


where referring to Figure 1, 

V*i is the displacement of the fuselage skin at ^ 

* 

Is the pressure on the inner surface of the skin at X 

is the displacement of the trim panel at x 

is the pressure on trim Inside surface at H - (cavity side). 

The coefficients are related as follows: ' 

^ll* t ) ^ii 8 a H * H 
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The representation above is a simple one that will lead to a 
tractable problem. It is also felt to be a reasonable one 
because of a number of factors. Its simplicity lies in the fact 
that the pressures and displacements occurring at each location 
are assumed to be related (across the trim) on a point-by-point 
basis. This may seem to be a very restrictive model, but since 
the transmission of acoustic waves is to be considered, the 
pressures and displacements will vary rather slowly with Tc . 

One could argue that for oblique incidence the pressure at a 
point on the incident wave front would be felt by the outer and 
inner surfaces of the insulation at different coordinate posi- 
tions, say 5C and X' as it passes through it and that the wave 
would travel a distance greater than the thickness of the 
insulation blanket as it did. However, the acoustic wave inside 
the blanket propagates in a direction that is much closer to the 
normal than does the wave outside the blanket because the speed 
of sound in the blanket is less than the speed of sound in air 
(Snell’s Law). Because of this the blanket does not appear to be 
much thicker for oblique incidence than for normal incidence [4]. 
Thus there is justification for the simple model postulated. It 
is important to note that the transfer matrix representation 
given above applies only across the trim insulation and lining. 
The fuselage skin (sidewall) is not involved. The sidewall 
response obeys a much more complex dynamical relation. 

The response of the sidewall (for harmonic excitation) is 

w, (JC > zz J ^(x')] d%’ } (1) 

where G(*|XJU>) is the structure’s Green’s function, the form of 

O «<• t 

which will be presented shortly. (K) is the exciting exterior 
pressure field. The integral is performed over the excited 
structural area. 
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The interior acoustic field acting on the trim lining is 


$(%'> = — Gl ? (%\x]tO)\Af z (X)6x } ( 2 ) 

where <*p(S\S';*>) is the Green's function for the cavity and ^ is 
the density in the cabin. 


The exterior field is ( exterior air density) 

P°C*') = (*) 4- p 0 ^J Op (*Y£;uJ) 

= + pr(*) i 


(3) 


is the exterior space Green's function, 
pr<*'> is the radiated pressure field on the exterior surface, 
Pb.K'l is the blocked pressure field on the exterior surface. 


The three equations (1),(2), and (3) and the trim transfer matrix 

(which consists of two equations) forms a system of five equations 

• 0 

in five unknowns: pj 1 , , y^, V/|> an< 3 |? r (or jp ). The blocked 

pressure field is determinable once the geometry of the structure 

is fixed. It is now necessary to solve the system of equations to 

determine the effects introduced by the presence of the trim. 

Before doing this, note that if the trim is removed from consi- 
• • 

deration, and W, = \W*. reducing Eqs. (1), (2), and (3) to 

the previous set solved in References [2] and [3]. 

To begin, pj - and v/^are to be eliminated as variables. Prom the 
transfer matrix 

J 3 ! °*l| 06 \ 

= + (a 2l 



Therefore 


f! - = +Cv.w, , 


(4) 


where 


r — . r — dr-*. # 

P I +«*<*„ ' w 

Also 

w i= "ixW 2 

giving 

... _ Wi auR 

w, — — — > 

or 

= c p h l + C w v/, , (5) 

where 

c p = ^ ; g*=-±- • 

An An 

At this point, it is emphasized that in Eq.(4) the coefficients 
and C^are subscripted and in Eq.(5), the coefficients cT and 
C w are superscripted. These coefficients are found to play 
prominent roles in describing the trim effects. 

Continuing now, Eqs.(^) and (5) are substituted into Eqs. (1) and 
(2), and after using the defining relations for the Dirac delta 
function, i.e., in the forms 

W t CZ) = f w,^') > 

and 
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it is found that Eqs. (1) and (2) become the following 

J £$(X-X) + c w w t Cx')dx' 

= JG (xlXjt*)) [Vtx> “ Cp dx'. 

J D(x-x') + pu ) a cf Qp(xlX)u))J j?t(x^dx 

— -puT J*Gj p (xixj uS) V^, dx . 

Eqs. (la), (2a), and (3) must be solved simultaneously for the 
modal displacements of the fuselage wall. 

Consider the left hand side of Eq.(la). Let 

*,<*) = £s't ,s <*'> > 

S < 

where £ s is the modal coordinate and Y Cx) the associated modal 
function of the fuselage wall. Then the left hand side is 

•+• c w G(S|x;w)J 

r — v « V — ' f \ (x)V'(X) . 5 

+ C wZ_/J 2-^ MvYrC^ £ s ^CX)cix', 

S 25 V- 

where the structure’s Green’s function has been introduced 


(la) 


(2a) 


Gi(xixjio)= ) r 


MrVr(W) 


with modal mass 


and receptance 


M ,=/ m 


- 10 - 



where is the resonance frequency, l| r the mode's loss factor, 
and m the surface mass per unit of area. 


Using the property of the orthogonality of the modes, the left 
hand side of Eq.(la) reduces to 

E t ' ♦ -ffcj • 

is a dimensionless quantity. 

Therefore Eq.(la) becomes 


E[» 


Cw 




;] *W,= f Q(x, *')*>■>[ *>V> - c f j>*(x')] dx'. 


The above is multiplied through by and integrated with 

respect to R to obtain 



= -^r [/ — c r J~ t> l ‘csjf«')dx'] 


But according to Eq.(3) 
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Therefore 



Now consider Eq. (2a). The right hand side becomes 


( 6 ) 


_^u) 1 c v 'J GJ p (Xix'}io) J2 ( s V'c® <bC. 

s 

Multiplying Eq.(2a) through by V(X) and integrating gives 

j [sc*.-*') +• j>ilx)d*dx' 

— — pitf* - c" jVcx) y q P (xix'ju»r. f s V'cx) dx Jx' 

= - r ^E 4 ffq r (K |X» JxJx' 

= -r l c-E rv. , 

s 


where being consistent with previous definitions [2,3] 



= QfCXfXju)) K)lls%X>tAx<Lx' . 


The left hand side can be rewritten to yield 

r P r +• c' r * fftfc ^p(xix-co) j»*ck-> dxdx' 



u>V" 


E r s e s • 


(7) 
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Return now to Eq. (6) and consider 



r£ = J f’isi) 

-J V f &>[p o^'J «**'] J* 

= C s f[ Qp(xH.'j*o) ¥(& V'(X)dx'Js , 

S 

or using the previous definition of J r *fco) as the integral 
appearing in the expression above [2,3] 



(8) 


Now substitute Eqs. (7) and (8) into Eq. (6) 


S 


-yh^’Ens,] 

w 


“ c r ~ fpfc, • 

Collecting terms in ^ gives 

{-**[»*!&] +^[p ,r + f c pC -I-]}^ 

^ 4 E[f'^f c r^I rs j4 

s*r 

+ c r c? r 'ff 'P$k’> (q p ft |X;<o) |> 2 l (X) J* dx ' 3 -fp ^- ( 9 ) 
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Consider the integral in Eq.(9): 

J j tyc*') Gp(xiKj«>) <Ax<Jx , 

with the cavity Green's function being given by 

g, = r . 

where 9 «<*>„ is the acoustic modal function for the nth mode of the 
cavity and is the complex eigenvalue. According to the 
previous normalization [2,3]: 


JjJti j7 = 


The integral therefore becomes 



jnj£ nWjhtf) r , _1 

v k*) L Z_, . 


where 


K ix) • 

m 

with 7T m = normal coordinate for interior pressure (dimensionally 
is equal to pressure). 

Now the surface integral 

J ‘Pn^ > 

is identically zero in almost all cases whenW^n. However, the 
surface area over which this integral is evaluated may not extend 
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completely around the periphery of the cavity since that area Is 
the transmitting surface. In such a case, for a few acoustic 
modes, the above integral will not be identically zero. However 
it will be small in magnitude relative to 



Therefore, the following approximation is valid 


’*»]**’ 

€ *a /C ( * >elR f 

V V (" X n- kX ) ^ ’ 


where 


and A is the transmitting area. 

Substituting into Eq.(9) gives 

p r] K 


E [«r- Vr 1 "]*. 


+ i—t L ro 

s^r 


‘ B - r * aa 
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Now reconsider Eq.(6) which can be put in a similar form 


+ jrs ^ 




(ID 


Note that in the above use has been made of the relation 

= / V'Cx) fi (X) lx = f tp- (X) vt>n (X) Jx 

n 

^»A-f (n,r) * 
n 

It is now necessary to eliminate the 77^ from Eqs.(10) and (11) 


Let 


- M,V r [l+-^J + c p c w ^r r ] 

a rs = w^j^ + CfC^^r] 

l _p .a. ■f'tn.r) 

b rn - Cpc^co- • 


Then Eq.(10) becomes 

°' rr ^r +13 A «f s +53 b-«7T n = — ^ , 


or 


s*r 


n 




(12) 
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Let 


c rr = - tAr \ [ \ +• J +• 


x *»- -T'^” S 

^rs — po*^ ^ 
d rtl = — CpAf'Ov) * 

Then Eq.(ll) becomes 


. c J9 + MM = tn.} • 


From (13) 
or 

(7r} = -[d]'[c]{f}-[d] ’(t ) * 

Substituting this into Eq.(12) gives 

[*] W - [»P !<#) - [► P Snj - t w • 

or 

{ f) =[[«]-[ b i4c]f[[ b ][d]:[r]]{irj- 

U4) 
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From this form, it is not obvious as to the influence of the trim 
on the sidewall response. Thus, the terms composing the 
coefficient matrices [a], [b], [c], and [d] must be examined. 

First consider things that are known. Since coupling of structural 
modes by the acoustics is negligible, generally 




J rr 

can be included in 

the loss factor >\ w . Thus it follows that 


V 

C 


rr 




and 


C m 2s o. 


Therefore [c] is a diagonal matrix (very close) = C^c-J. 

Also X and X are so small that [ a ] is diagonal: r*a.J 


-rr 


-MrX[l + -£%-] 


a«ss o. 


It follows that faJS (very close). 


Then 
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But the inverse of the product of matrices is determined from the 
relation 


(ABC) - ' = cW. 

Therefore _j . _j . 

= l> Jl-rk. } = Wfrj- 

This form is quite similar to the case where trim is excluded. 


Without trim 


With trim 



where 


The solution is 



, I 

£ t + m YrJ 

— Qta. Aft 


^1 X!/ ’ 

giving ^ 

r t 

which since ^*0 for t+r reduces to 

— -~Y2 ^<*> *»■•- • 

i- 
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(16) 


( 17 ) 



To 


get 


from 




n 

Eqs. (12) and (13) have to be reconsidered. Either can be used to 
obtain the 77^ . Prom Eq. (12) 


or 


H{fM b ]W=K,} 

HM+MW=fr n ,} 


But by Eq. (15) 


W-NfrJ 


that Is, 


where 


W=NKU + ( f <} 


1^*1^ I otTr ^i I ■ 

This gives 

H[HkMO] + [ b ]M 




or 




W— [*-t«J(0' 


1 .e» 



Prom Eq . (5) , 




W z (x) = C^(K) -f c* w,(x>. 

Thus 

W 2 (*)= c f { <f> ]{ tt) + c w { ^ }{ £ } 

=-^4 b l*j{o + c"w T {f} 

— C v 'v/,(5C>, 

in the limit as acoustical coupling of structural modes is 
neglected, i.e.. 

The interpretation is this: while it is necessary to have a 

« 

precise calculation of 77 *^ for purposes of determining 
the values of 7T n have little bearing on the calculation 
of and for all practical purposes can be set equal to zero 

for the W^CX 1 ) calculation. 

This now allows a calculation of the power flow into the cabin 
using 



NAf t (x> = C W W, tx) • 

Power Radiated into Cabin 


( 18 ) 


The cross power spectral density of pressure ^ at position X 
with trim lining velocity at is (for a one-sided spectrum) 


j(^) P*v ( * ;u,J • 


- 21 - 



^ and are Fourier transforms of the truncated random 

quantities (* = conjugate). 

Now 

ViT U;“j)= i.toW^.(K;w) . 

Here the definition of the Fourier transform is 

3K )=/" 

— OO 

Therefore 


Spiv^Kl^iO) = 


lim .2 <-<0 
T+o* T ^ 



But according to Eq.(2), 


Thus 


Pi T (*;'°) = “ p* 0 ”/ I*;" 5 ) 


= -Lpua 3 


T->~ T ZTT 


X J Qp(xlx'] w) W^x’.'u)) dx' 

= ~'^f ^VV 2 ^>)W*sM]dx" 


= -L p to 3 / (q p (x |X j ui) S^(x"ix 1*J) dx'.' 


(19) 
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Now let X -►* • spectral density of power radiated internally 

is obtained by integrating over the transmitting 

area. 

mt f 

^t'AA ^ “ J dx . 

This quantity is complex. The real part is the real power (i.e., 
cospectrum of power). 

Therefore 


VVr*j M = —A j Qp(x|X*-u)) S Wj> (x'l*jto) dx'dx . 

n' 


But according to the previous finding 


W : t K,w) a c w W, t (k,ujV 


Therefore 


= -p_»oo f (ilf) ^2. T W ZT (XjU>) 


it- ^ fe) cXCK«) 


nV/ 


SyyjCx^lXjtO) , ( 20 ) 


- 23 - 



and 


Wf^Cu)) * - jc w | ^ p (xlx , )u))S W| (xix;^)clKJx , . 

KK' 1 ' 

But 

W, T (K,^>) = -J2 V Of) <*.-,■ fpjjcio) . 

So 

^W^x'lXjiO) V'CX)^*) 

r /* 

x S^(xiKju)) V'(k)^cx , )4xdK , > 

giving, upon introduction of the joint and cross acceptances, 

m = - A f co 3 | c w | 1 ^ * rr «u[ff G p (m';u>) - 

r > 1 jtx' 

x V /l Cx 1 ) V \%) cis'd*J j^e (**>) to) A 2 

<A ptO C A (w) • 

l r,A 

Neglecting cross acceptances, this reduces to — 

pui 3 C*| A i Sp b| (u5^|q^.| I N) j^u>) 

r 

Now 

i (w) = — — r^Cim] + ictucni)] 
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Therefore 


*.[■ 


w:>] - f^A 1 ] <-Pb^e 


M r (Yr+-C w /m)| 


= r *AV-i.|c-|‘EfE#%.fi, 0 


Jr< w > 




_ 2pcoM 4 I W |*p «, X„A„y^ ffi r) j r \«M c (u5 v 

It follows that the band-limited power flowing to the acoustic 
space is 


Aio 

s 
v 


• 2 


_ P C 0 A 4 |-,-W | l /•* , ^ JlHtlO) 

fc I S FmHL -nr 


(21) 


- M? 

x n €n *l»i*f (Wif) r to*du> 

V w " -L J ( 1 “ Yr+ ’ 

where use has been made of the relations 

A„ = /c 0 > A. n = co h yj n /ic„ , 

X „ — A „ — * A- n • 


and 
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Let 


Then 


and 


Cyy ® Cyy + A Cyy f 


Y.-h% 7= 

w\ ^ 




Y,+- 
r m 


= u) 


I C-w I 2 c5. x. 14 

l| r u)r +V| r tt) r 4 


+ 


2 r 

w _ 2c w 


Since in general 




I ^ w | 


m 


the above reduces to 


Y r + 


C w 

m 


• [L *V\J ^Lm\d« 


2.Cw Hr , 
— + 


m aJ- 


Now set 


W x — Jfnl 1 — ^£l Jk. 

> “ m'u)* m * 


This leads to the result 


Yr+^f- 


m 


3 

v» r - i y 
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and Eq. (21) reduces identically to the form of Eq.(3), Section 
(1) of Reference [3] except for the presence of | C** I* • It 
follows that 


w.-_ = 


<f>tl> 2irA’ 




fV 


Ifl 


f 2 Ui’utf r [ C r-Cn \ 

Mr U D„r U 4 ]' n 


+ f ZCn(b f -b B )- bn(Cr-Cn) \ ftr& 4 - an -(. / C n -C. r_) | n 


2L c r ~ br (Cn-Cr) 


where for ft or f | 




-j arc+an h J 


(22) 


i . f IO^C w /2) 4 U3 4 4»bj(l^C w /2fa3 X 4-Cj 

• nj iM I. a A L , .z . ^ 


with 


|(l- C«/2)V+ bj ( I- +- Cj 

_ .x _ j_ ( 2 +Ci U )\o’’- 4 tOj i -i ( 2 -C lo ')\o ,: - 4 u>J 

arcxan •. = tan =• — tan ^ > 

J 4>|j uJj 4-rpujJ 


Dnr= (t r -c„f + (b B -b r )(b n C r -b r c n ) , 
b n = — 2 -toJ ; b r ==- 2 «>r > 

C-n-ujJi i Cr- i»>r(i+-^r*') • 


The above is identical to Eq.(5) of the Reference [3] except that 
J | Z multiplies the result and Yj{, replaces YJ^ . 

Now closer consideration needs to be given toC^and C W . 
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First consider the multiplying factor |c W J* • Suppose the trim 
Is a simple model of Insulation plus limp mass as shown in 
Figure 1. The transfer matrix for the insulation is [refer to 
transfer matrix derivation in Appendix A] 
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When coupled together, the transfer matrix across the trim is 
found using 


This gives 


h . _ h 

Vl J Vi 

t«»vUKon wl* 
output Vtuss 


ft 


\ « 


- 



r n 


_ b 


1 

tu)rw t 


% 

1 

O 

ft 

V, 

✓ > 

1 • 

0 

» 


_£ c 
L w C J 

V, 

> 


(23) 


Also 


output 


V, = -au> W, • 


ink 

intuUtion 


The transfer matrix of concern here must be of the form including 
displacement. Thus 


f ! 


Vi* 






P 


I 

0 l 


C LwWS 

-LI c 

J Ltow ^ _ 


ft 


(iu)Vfs-no*m t C) 

-ii. c 

6t>V^ 


ft 


(24) 
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Now according to the definition 


or 




w 

CW - S 


C,S and w are complex. This leads to 

i^i 2 W w* 


CW — LtonttS Cf Vs/*+lu)mtS* 


Iwl 

•z 

icriw 

a *f toViJ < 
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A Preliminary Comparison 

At this point it is useful to make a comparison against some 
measured data. Consider a trim such as that of Figure 2 similar 
to that used in the 7^7 aircraft. In this case, the trim panel 
mass is “ 1.76 kg/m^. It covers a 0.102m (4 in.) thick layer 
of Owens-Corning PF-105 fiberglass weighing 9*61 kg/m3 (0.6 
lb/ft3). 

From Appendix A 

Y rr oc — l k , 

where 

K — Z-'R £ . 
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is the acoustic wavenumber in the material. The attenuation 
constant is 

= oc /Z.(o<\ , 

r ' 

r 

and the wave impedance is 

W ’ = ~ "Zs ( ^/ yy ^‘" s ' C€ 0 ' 

At high frequencies, 

I Vs / 1 — > 42.0 k^/^-scc . 

Also 

V/ = |w|=^ = |Z s |c" i6 - Z: s 

The properties of the fiberglass are given in Figure 3. 

The prediction for Jc^i is shown in Figure 4. On the basis of 
these results the trim panel vibration isolation needs to be 
examined at high frequencies, otherwise, a reasonable prediction 
is afforded with the model. 

Other trim configurations can be easily generated. Suppose the 
insulation is replaced by an airgap. Then 

W =. p e c. 0 

} w/c D . 

Cppo) 

C | 2 +- w a m£ | S x 
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or 



1 -W 1* - 

1 

|C | — 

1 +■ (*mw % /ac 0 ) 1 

This is simply the theoretical 

"weight (mass) law" 

spoken of in Beranek and Work 

• 

1 — 1 
in 

i i 

At certain high frequencies 

9 

go* ZkL.-*'© y 

si«2.kL -*■ 1 


1 

1 4m 



J “\f A,/ 

Since Wm t / p,C e » 1 , 



which is mass law. 



This also looks good. 


attenuation 


As has been shown then there is appearing in the theoretical 
results, a trim coefficient which will be called y T^~lc VV |% 
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This is the "add-on" attenuation of the lining and insulation. It 
should be noted, however, that the multiplier n is not the only 
effect of the trim. There also appears added damping due to 
, as previously discussed. 


There will be occasion to consider the lining dissipative (there 
being some flexure of the lining but with the lining stiffness 
small enough to ignore). Then the trim transfer matrix takes the 
form 




(Lu)tti t — 






> 


where 



is the trim loss factor. 


In terms of displacement, which is required in the present formu- 
lation, this becomes 





✓ 



(coVi^-fL. 


L° 


1 

- 

- 

N 

A»i 

4-IX 

P- 

*X| 


w, 

✓ 




C UWS 

-ii c 

L.W\N K 


\ 



V/, 


✓ 
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Therefore, for this case: 

°u = C. + tonuv| t 5. _uom t i , 

a (1 = io*m t C 4- i_ ( ui\NS + u > 1 C.) > 

„ - _kJL 

“•*! - uAN ’ 

— C . 

The trim coefficient is always determined from the coefficient &|| 
C w is determined from the remaining three coefficients. 


Band Limited Power Absorbed on Cabin Wall (with trim) 


This result is taken directly from the previously cited reference 
C 33 - The required expression is from the appendix of that paper, 
Eq.(A5). First the conductance, £ , looking into the trim is 
considered to be independent of X (do not confuse the conductance 
£ with the modal coordinate of the structure, i.e., ). 

~ dv ) * (26) 
< > 5 ^ is the space-time mean square modal pressure. 

v W 5 7lrID < « < i • 

' n * 

when the normalization of the is taken into account. 


where 

Thus 


-3*1- 



In the case where 


i 


Is a function of X , 




e n< 



( 27 ) 


where 


■/, 


Is always evaluated over the absorbing surface area. 


The conductance can also be related to the imaginary part of the 
eigenvalue of the acoustic mode, or, alternately to the loss 
factor of the acoustic mode. From [3], if £ is independent of X 


■A-n 3 * 


k 

Z\n 



t 


If not. 

In terms of loss factor 






Thus substituting Eq.(28) into (>27) gives 


( 28 ) 
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In this report W^ n is always used in the form of Eq.(22) and 
^diss in the form of (29). The relation needed to deter- 
mine the effects of the trim in terms of the absorption from the 
space will be Eq.(28). In other words, the acoustic loss factor 
of Eq.(28) will be calculated once is determined. Alternate- 

ly the loss factor can be measured. For the bare cylinder case 
(or a bare fuselage) measurements will usually be required. But 
for cases where trim is present the analytical determination of 
is desired. Therefore, the calculation of the conductance £ of 
the sidewall is needed. 

By definition, the conductance is the real part of the admittance (3 


i- Re(£)= R«(£-Jl*) 



= , (30) 


where V and p are the trim velocity and the pressure on the inside 
(interior side) of the trim. 

This can be written 


13 - - '(.tOpe-aW/p , 

where W is the trim displacement. 


Suppose 

v = Ve- l,ot 
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Then v/ P = v/Pe- 1 -* where V and P are real amplitudes. If 
<P=°,0= real s=. f = ^.(V/P)- When p is maximum V" is maximum in 
the direction of power flow. Be definition, is positive 
inward toward the interior (as is ) - refer to Figure 1. 

However the conductance desired is that looking into the panel. 
Therefore V=s — V x (or W»— Wj, ) . Thus P is determined with the 
relation (using r-t* > 


Now 


and 


0 - Ltopc.v/j/pi . 

% 

-t-a.j.v/, , 


(31) 


w L = a t| Y, +- Ail'll • 

Define the impedance looking into the structure at )C as 

p\ = V = ± 

V, “lulW| tu>W| 


(32) 


Again the sign change is necessary to define an outward velocity. 


Then 


0 — tiopc, 


&XI AtrW t 

L A «t ft +• 


= P c ‘ 


-•to a x , z, -t- 1 io a zx 


iu) a 


n 


L, -b CL 


IX 


( 33 ) 


This is the general result desired. 
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Before continuing, a check for correctness is needed. Let the 
impedance z' be calculated, where 



iu) flu 


For the initial trim model it was found that 


a„ = c 

q 12> =s LuiVs/S +■ 

a = -UL 

2 ' ujVN 

0.2.2. ^ * 

Therefore 


/ __ ZidcoC^^mtS/w) •+ (lajyys -4- 

( La) S/W) -4- cu) C 

Dividing through by La) (numerator and denominator) gives 


Zt (C-i-qjnVfcS/w) + (WS -tom-t C) 
S,(5/w) + C 


W 


*2, ( c- Lu)m t S/w) +(VMS -*cu)iv\ t C) 
Z.,S H-NVC 
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Prom this, it is seen that the impedance looking into the trim is 
of the correct form. For instance, if -+• O f 

z'-\N 

which is the correct form for the insulation alone [6]. 

Now return to Eq.(33), repeated below 


Cz, -t-WS 

Sit) +• V/C 


/?= 


r c • 


— oo a.- u «<- c to CLx x 


(33) 


u to a u +• Ol\z. 

2.,= j^ytuivv, = impedance looking into structure. 
In general local, and therefore 


d (x) = pc. 


— ‘tO (x) 4“ LcO a XX. 

iWOLyy^ C*) +* 


Now 


But since Q is sought, one can for the present consider 

f>° a) = f> r ts.') , 


i.e., only the radiated pressure on the backside (outside) of the 
skin is of concern. 
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Recalling that 


{V*) = p<£ L J Gip (Six') w,cx) dx, 

and noting that the radiated pressure effect can be included in G 
as increased damping since the reactance of the air is negligible 
compared to that of the structure, for the present purposes 


\N t t>U = —J" 6 (xixju>) p, c lx5 d* 1 


M s Y s cw) 





V'fc) 

Ms Y s 




This leads to an approximation 





s 
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One model is obtained by letting the pressure be uniform 
( constant) 


» • 

ft w - \>\l+S * 


leading to 


L - MsYsM - 1 


(3^) 


it L 

Since for a given interior mode, pj follows ^ , it is also 
possible to write 


% 

1^*)= C„ <£>„(*) 


> 


where C n is a constant. For such an approximation, another 
structure impedance model is found: 


^,(>0 = 






( 35 ) 


and the backing impedance depends not only on the structure, but 
on the prominent responding acoustic mode as well. Thus a precise 
calculation of is a complex matter. 
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As will be seen, if the trim is fairly limp and absorptive, 2, 
need not be very accurately determined to obtain Q. In fact, in 
the present work, a simple model 


2 , — |OC d - LtOvv\ 


has been used. 

Noise Reduction Calculation 


This is obtained by setting Wi n =Wi',„ from Eqs. (22) and 
(29), solving for the mean square modal pressures <; 
one-by-one, and adding according to 


H < P" 


Sit 


The result is found to be 


< _ 8 A 1 fm A_\ z a | »«■ 
<pi>v “ c » w v i l4j '»i c i c ’ 

C J r (u)) T(«l, r ) 

Ml D, 


co : 


+ 


n r ,w ‘ir 
ZCin Cby — b n (Cr~ CiQ 


Cr-C v 

4 

otrc-tan, 




(36) 


where the space average mean square pressure •<, is (for 

the reverberant field) related to the mean square blocked pressure 
by 


< p«> s ,t =<)>:,>/ z - 


-1*2- 



r i= 


= trim addition transmission coefficient. 


= modified structural loss factor. 


V 

The noise reduction is 


NR= -lolog ^^^ • 


High Frequencies 


Without trim, from the Phase I work 


W. rr 
v, «n 


At. 

I 4f<= 


+ - F 


^■etz^Lo'H 





Hirr 

IT m* c 0 Z_fL J 1 *'' 0 ) 

r<.A«o J -A 
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In the above expression, the second term represents a power flow 
through the panel (structure) outwardly. The transmitting area 

is A. 


In the case of an added trim, there are three areas of concern: 

1) the transmitting area with trim, which will be called A*. > 

2) the absorbing surface area S, and 3) the transmitting area 
without trim A x . The total transmitting area is 


A = A t + A' , 


and the inflowing power takes the form 





A'lg 8 A A' 


/*- i / 


TTlrv\ 


- cbE [W 


la -1 


< I* >s,. 


f-iAto 


.% a .% + jaa; r ^ > . 

n-(T°) + »«*■ ^L J J j =a 


(38) 


Tj^. is the sum of the average external radiation loss factor, 
i.e.,Tj^^ , and the average structural loss factor "f^ is 

the trim transmission coefficient, and / 1^ is the resonance 
transmission coefficient for the diffuse field case. 


R. = ^i/p^r) ■ ( 4 P c «/ A ) 


In this expression is the modal density of the fuselage 
structure (modes/rad/sec). 





The radiation loss factors for a diffuse field are given by 



2 P„u>A 

1 TlMCo 


< j 



2p A 

TWCo 




> 


A. 

where <J r-'V is the joint acceptance averaged over the structural 
modes resonant in the band. 


Let = field incidence transmission coefficient for mass 
controlled panels 



32 * A 

IT 1 m‘ 


E 

r<Aui 


as defined in Phase 



Then Eq. (38) becomes 


<&>*.* 

= 4^ V^) A t + (T R +T f )A'] < ^> st 

4 f>C„ [ IjgJ t. + ^]< ft >s,f (39) 

1 Eq.(l6), Ref.[l], p. 20 is incorrectly typed. A should not 
be squared. 



The power absorbed is 



( 40 ) 


S includes the area A covered with trim and any other absorbing 
surface area. Equating Win to \N ii%s gives the desired high 
frequency result 


<le>s,t _ cxS ■+ + Tr( >)r/V)“j )] A' 

< P l ‘ >s * % ( T+ -i-Tr) At + ( )A' 

Letting gives 

< +[l> + Tk( V\r/ r l^d)]A / 

< f^s,* TA' 


This is the fundamental result with trim present. Note that if 
all transmitting surface is covered with trim. A* = 0, and 

< pe >g .t _ 5-S £;S 

<ff>s,t T t TAt “ tfct-A 

On the other hand, if trim does not exist, = 0, and 

«S + [t-f + T r ( y\ r / A 

< ft> s ,± r A 
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which is the Phase I result without trim 


Finally, oC is estimated with the relation 

s > ‘ <1,2: 

If is calculated (say with trim only absorbing), Eq.(28) is 
applicable. 

3.2 Influence of the Internal Radiation Damping 

It has been found previously that, in the relation for the 
band-limited power flow to the cavity, there appears an integral 


f _ 

co’d 

u> 


J 

4to 


z 

M r Y r -pu>*-( 



In Ref.[l], this was evaluated after redefining to be the sum 
of the structure's dissipative and external radiation loss factors 
and by neglecting the internal radiation losses. This allowed the 
term in the integral with J^to be incorporated into Yr and the 
term with X to be neglected. When this was done, the above was 
reduced to 

f 

jjXn-kTM?|Y r | 4 ’ 

where 

I a* - kT= ^ [( , 

and 
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with 



_ _ 5+r*c+\f ft. 

V 


r*.d 


The integral becomes 


M r l -;< -L {[| - -£r] z + 1«}{ [( - -^fJV 1?} 

which appears in Eq.(3) of Reference [ 3 ]. It has just been shown 
in Section 3.1 that when trim is present the loss factor r^ r above 
is replaced by Vjr where 


/2__. |Cwl 


m 


to*- 


+■ 


Now for increased precision, the question must be asked: what if 

I rr 

(UJ) is not neglected? In this case, 

consider 


I kx - p^rf 

= K | x.u» - ru 


where 
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Consider a case where u) n £:^(*) and close to u) r ; then consider a 
single acoustic mode, say mode n . 



P'A* T rr 
1 



p«O l A 4 €f.J f Cn 4 n 0 uJ* 

.)+ ‘-'In 

1 

M r V*o£ 

[(’ «•*) 

2. 2 1 



This becomes 



^>a)M Z € n cJ f(n,r) 

[-& 


M r V 




e n e£ fcn,o 



Now it can be shown that the non-dimensional form 


y _ P * 0 ^ cj <c - 
7 M r V 


I 


For instance, for the stiffened 0.020” cylinder of this study, for 
the two modes 


n =( 1 ,2,l) ; r =(2,Z) in 630 Ha band 
f'cn,r)t= O.IOZS 
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and 


CU n — 67Z X 2-TT ) Ud r = 660 X 2.TT 
p = 1.2.0+ k^/wn 3 

00= 630 x -2.TT± Auj/2. ; AtO= C^OO 
M r = mA/4 ; m- 4.21 kg/w’- 
c.„ — 343 rti/ s e c. 

A “ k 8g w ( Sufldct cyliwdtp) 

V= O.Z+3 iTi 3 
y= o.oig . 


As 00 ranges across the band, the reactance term 
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ri-^i 

L ^nJ 


l , ~w) + 1 *. 


takes the following values: 

it 1Z 


If 


»»<"V ['-■$] ^ 
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and 


If 


r 

K| * and again 
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[l — 

If 03 =*£■ UJ n (very close) 

^ ~ ^=>0. 
rjn 
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Thus the reactance term is always small (in air). This gives the 
approximate result as 



Now the structural mode appears to be damped more heavily, by an 
amount Kjy_ given by 



( 43 ) 


For any given 
since and 


spectral component at frequency u> 




are small 



pud*- 

Tv 





we see that 


is dominated by the reactance term i -c , unless tt) is very 
close to u) r . When a) is close to uV , the resistance term 

becomes important. Moreover, if simultaneously u) M lies close by, 

ft n 

r|f. can be significantly large. Thus, a frequency independent 

form of the resistance can be used and the following approximation 
is assumed to apply 
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- ^)- L (nr+V')] 



where 



When trim is present 

r \*-*- y \r • 

3.3 Transmission of a Tone 


In the case where trim is present, assuming hereafter, that a 
floor will always be considered, the tone transmission calculation 
of Ref.[l] can be generalized to obtain the solution for the 
inflowing power to the cabin, i.e.. 


v f V 


n n tr Lj us 


11 1^74 // C P (xu) i'm'i'w dxdx' 

y r 5CX 


3 


[( 1 - %h *ih] [( 1 - §)+(*]( + 


(44) 


n 
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y|£ and are used to include the trim. The above result is 
consistent with Eq.(33)of Ref.[l]. To review U) e is the frequency 
of the tone. In the case of propeller noise, this equation will 
hold for the blade passage frequency and also for each of its 
harmonics. The correlation function describes the 

blocked pressure field over the fuselage (cylinder) and is defined 
using the following relations: 




Cp b i(xix;u))== 2 Cp^(xlx;u)) ) o±c , (46) 


C F(> ^I<U))= (-T) 

where is the average cross correlation of the 

pressure over the blocked (immobile) fuselage. 


x 

=^ - 4 r [ T ^ Cpbi^S^] Re [ j\,C>s -t + n] dt- 

** (48) 

A major task for Phase III is to develop a suitable representation 
of for Propeller excitation. 

The power inflow to the cabin wall from the interior is as before 



( 49 ) 
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Equating Eqs. (44) and (49) gives the interior space average mean 
square pressure attributable to the nth acoustic mode. 

The interior space-average mean square pressure is given (again) 
as the sum of the space average mean square modal pressures: 

< f?V =E < £>*,*= • <50) 

n 

The average interior pressure ultimately must be related to a 
significant exterior level, say the maximum occurring nearest the 
propeller plane. 

Note, since 


< 



/<£» > 


where <. |?*>Is the mean square modal amplitude, the modal ampli- 
tude is determinable from 


<f£> = €„<^> st , 

and theoretically, the mean square pressure at every interior 
point can be obtained with 


<K^“ ) > t =E <5l) 

n 

It is emphasized here that a good point-by-point prediction may 
require better input data than can ever be generated, but the 
point-by-point prediction might be quite informative 
nevertheless . 
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3.4 Modal Properties of the Cabin Space (Cylinder with Floor ) 


In Ref.[l], the acoustic modal characteristics of an 'ideal 
cylindrical cavity were presented and various quantities were 
calculated with the modal data for use in making the interior 
sound level predictions. The modal data consisted of the mode 
shapes £ being an interior point, resonance frequencies (J n , 

and the mode normalization given by the constant In the 

present case, consideration must be given to the determination of 
the same quantities for the case of a cylinder with a floor parti- 
tion defined by the angle Q 0 as shown in Figure 5. 

In the ideal cylinder configuration the modal properties can be 
determined in closed form by an analytical solution of the wave 
equation, subject to the appropriate boundary conditions. This is 
possible because the wave equation is separable in cylindrical 
coordinates and the boundary conditions can be expressed in these 
coordinates. When the floor is present, it is no longer possible 
to derive the mode shapes analytically since the boundary condi- 
tions are irregular. Thus it is necessry to resort to numerical 
methods. There are two possibilities, either finite differences 
or finite elements. Since in the present case, the modal charac- 
teristics in the axial direction are known, a two-dimensional 
problem remains, and thus the finite difference technique, which 
is the simpler of the two is chosen. In two dimensions, the 
finite difference calculations are fairly efficient (fast). 

To begin, the two-dimensional problem is solved. Then the axial 
modal information is factored in. Next the normalization of the 
data is defined. Finally, use of the data in the calculation of 
the various quantities required by Eq.(36) is considered. 
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Finite Difference in Two Dimensions 


In the cavity (cabin), the Helmholtz equation applies. In the 
two-dimensional problem, this is 


where 


VV> + k> = O, 



The central differences for the grid shown in Figure 6 are: 




Rtwt,n Pw-i,n 

ax ~ 

2. AX 


Rw/t-n ' Frn, n-i 

ar 

ZAY 


Pw+I,** 2 Pm,n "t" 

ax* 

(AX)* 


Rn,«+-i r“ 2. P*n,n "t* Pm, n— | 

dY l 

(AY)* 

and (f) are 

taken to be synonymous. Let 

ax = ay = h. 


Then the Helmholtz equation becomes 



Eq. (52) applies fop all points defined as interior (I) or boundary 
(B). 

For boundary point B, exterior points E are required to write the 
finite difference equation. The total number of unknowns is the 
sum of all I + B + E points (Figure 6b) . The number of boundary 
conditions that must be used is equal to the number necessary to 
eliminate all E's. The resulting matrix for determining the 
eigenvalues and vectors is 

[Alfl-Fffl -*{$}• 

where 

V = k"h l . 

Given the eigenvalue X , the associated resonance frequency is 
obtained with 

t Q'W . 

Ci ~ K ' 

or 

and thus 

x _ c«rV Xt 
* L ~ ztrh 

^X^is ordered such that X» u _, ^ X^ <, X.^_^ ♦ 
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This is the resonance frequency for % =0 -% = axial index, 
i.e. , 

kj= ‘fcii'/i-. 

The resonance frequencies for the cavity are obtained with 


leading to 


k 5- k i + 1- 

% = “iff [ k v + T^] 2 • 


Boundary Conditions 


The boundary condition is that the outward normal gradient is zero 
(the wall admittance {3 is assumed to be sufficiently small to 
allow this assumption). Thus, referring to Figure 7 


d±) 

dr /9,<e< itr-e. 


where 


dr 


This gives 


or 



dx ? 


£x 


cosG — 



sm0. 
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Now m,n are the X,Y indices: increasing m implies Increasing X ; 
increasing n implies increasing Y. From Figure 7 


a <f> _ a<ft . 14> b(j> 

a* " 3Y ' ^ "aX 


This gives the boundary condition as 


i4> 

ay 


cos 0 — 


M 

ax 


s ‘m 0 




where 


av 


Pm,n»l ~~ Pwill-I 

2h 


ax 


Pwvh, n P. 


**-»» V» 


2h 


Therefore, the boundary condition in the region 6* # < & < Zlt— is 


005 S,n ®>wn* (53) 


In the region — » 0 e <*- 0 <- , the boundary condition above can be 

applied by simply setting 0 whenever ^ * 
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Details of Calculation 


The origin of the rectangular coordinate system Is at the floor 
centerline (capital X and Y in Figure 7) and the cavity is symme- 
tric about X = 0, so only half of the cylinder is considered. The 
approximate grid spacing chosen is a/7. The exact grid spacing 
is 


h—Ga,— 


in’bwjr*! V»r+ £ *7 ( \ 






This gives an integral number of grid points between the floor and 
top of the cylinder. The location of the grid points is (NX>NY) 
at spacing h for 

NX= -MJ to +-NX > 

NY = 0 to NK , 

NK = Jj7 ( l +- CoS > 

Ihfecjr^j Part l 

The centerline of the fuselage has the coordinates ^ 0 } ACDS© 0 ^. 

To select the grid points representing the curved surface of the 
cylinder, for a given value of y-ny-$«, the 2 grid points are 
identified which lie Just inside the cylinder and Just outside the 
cylinder in the following way: 



and 


where 


and 
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The intersection of the circle with the line Y = NY-fia is given 
by 

X b = cl £ \ - (tose 0 — HY-6') 1 '! 

The last point inside the cylinder is (NXl,NY) , where 

NXI = p»<-t -(cosS e -NY G) i y /: J- 

The first point outside the cylinder is (wx+i.ny) . The 
distance of point (n^ny) from the boundary of the cylinder is 

aj l-^(cose 6 -NY-G) x + (N)C-O x y /Z | ' 

A test is performed on the 2 points (n*I,NY) anc ^ (NXX+ljHY) to 
see which point is closest to the cylinder boundary. If the 
outside point is closest, then it must be included as part of the 
cylinder. 

This procedure is repeated for successive values of MX for NY=0 

to hy**nk . 

This will not identify all the boundary points and adjacent 
exterior points, and a similar procedure is performed for 
successive values of NX* 

The line X«NX-Gq can intersect the cylinder boundary in 2 
places, depending upon whether 

NY-G< COS&o or NYG> cose«. 


If 

NY-0 < cos e 6 , 
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at the floor for 


NX • <a < sin ©o y 

the boundary points are given by NY — O • If 

NY- <k < COS e o > 

and 

NX-G. > Sin 6 0 , 

the Intersection of the circle with X=NX-Ga is given by 

fl-{cosB 6 - [ ! - • 

The first point outside the cylinder is (NX,NYr -0 , where 

NYU — \ — Par^ 

The first point inside the cylinder is (nx,nyx)- 

Again the distance from these 2 points to the cylinder boundary i 
compared, and the point closest to the boundary is included as 
part of the cylinder. 

Now in the case 

NY-G > cos e 0 , 

the intersection of the circle with X= NX-Ga. is 

\ ~ Cos © 6 -h £ | - (NX-G) 1 ] } • 

The last point inside the cylinder is (nx,nyt) where 

NYI = Xn+e^l Rvff {■^*^ c os© 6 -f(l-(NX'Q) i ) / ^J j * 

The first point outside the cylinder is (NX,NYX-H) . 

Again the distance to the cylinder boundary is used to select the 
closest point.. 
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In this manner*, the following grid points are therefore 
identified: - 

1) Cylinder interior and boundary points (total = ) 

2) Cylinder boundary points (total * Kl^) 

3) External points adjacent to boundary points 

(total = Hg.= n b + 2 ). 

These represent only half the cylinder, since the cavity is 
symmetric about X = 0. Therefore the total number of cylinder 
points is 


2 * K -(nk+0 > 

where 

NK - Iw+«<jra,l £ 7 ( 1-V . 

The total number of boundary points is 

2 H b •— 2. • 


Let 



= pressure at cylinder points 


ttiXl 


(interior + boundary) 


w. 


= pressure at exterior points, 
nexi points only. 


adjacent to cylinder 
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Using the recurrence relationship for all points I , i.e., Eq.(52), 
gives 

■ [ri]{Px} + [re]{p £ }= x{Px) • 

<n c* n 0 CWfXO (nixO 

For symmetric modes, symmetric about X = 0, when m = 0, 



For antisymmetric modes, antisymmetric about X = 0 



The matrices [RI] and m 
antisymmetric modes. 


will be different for symmetric and 


The number of points forming the boundary of the cylinder is equal 
to the number of exterior points adjacent to the cylinder minus 2, 
i.e. , 


n e = ■+• 2. . 


For each point 6n,n) on the boundary of the cylinder the boundary 
condition for 0 Wtt > 0 O is 




vn,n4-i 



and for < ©* , the above with 
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These boundary conditions give "b equations, and 2 more equations 
are required to solve for the exterior pressures iPe}- 

Two boundary points are selected, close to and 8*-3ir/4 , 

such that each of these boundary points is adjacent to exterior 

points in both the X - direction and the Y - direction. 

* 

At both of these boundary points , it is assumed that 


Rvi-i , r\ “ O* 


i.e., the gradient in the X - direction is zero (the gradient 
normal to the surface has already been equated to zero) . 


This gives a total of — equations, giving 

[bi]{p.} + [be]{6} = 0. 

(n e xwO ( M i*0 (*i e x (*e*0 

The matrix [BE] is non-singular and 

-I 


fe} =-[beJ[bi]{p i } 


(55) 


(56) 


Substituting this gives 


-I 


[ri:]-[re][be][bx] 

*0 (*L*ne) x v\ t ) (ricx tv *0 J 


M = X{Pt} ' 

(h- l xO 


(57) 


The eigenvalues and eigenvectors are calculated for symmetric and 
antisymmetric modes, separately. 
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The nodes. Hi, , represent only half the cylinder (i.e., X > 0) and 
the generalized mass must be calculated for the whole cylinder. 

The generalized mass for mode L Is 

V\i 

x CO) > 

j=' 

where h =Ga Is the distance between grid points. 

Is a factor, either 1 or 2 « 

1 for all locations on the centerline, l.e., X = 0 

1 for all locations on the cylinder boundary 

2 for all other locations. 

The boundary values of the eigenvectors are extracted from the 

L eigenvectors, and are defined for X? 0 only. The boundary 
points must be consecutive along the circumference of the cylinder 
from the floor centerline (0,0) to the top centerline (0,NK). 

The boundary eigenvectors will differ for symmetric and anti- 
symmetric modes for X < 0. For symmetric modes 

<£(x,y)= <£c-x, y). 

For antisymmetric modes 

<£(x, y)=- y)- 

Each mode must therefore be identified as symmetric or antisymme- 
tric when the boundary eigenvectors are used. 


C(j) 

r 

C(j)= 
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Apart from the boundary eigenvectors, there Is no need to distin- 
guish between symmetric and antisymmetric modes in the cylinder 
response program. The modes are therefore combined, and ranked in 
ascending order of frequency and the first 20 modes only are 
used. 

The first symmetric mode is a translation mode, with zero 
frequency and corresponds to the cylinder (0,0,0 mode. 

The eigenvalues, AJ, are used to calculate the resonance frequen- 
cies for the (three dimensional) cabin with 



where h is the grid spacing = Ga, L and a are the cylinder length 
and radius, respectively. 

The frequencies output by the program are 



> 


i.e., the frequency for o modes and(X=|in, i.e., for unit 
radius. The generalized mass is also output for a cylinder of 
unit radius, i.e.. 



G 2 o? 


y 


where O. =■ I m . 
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Normalization 


The acoustic modes $>n for the three dimensional cabin are 

where is the eigenvector for the two-dimensional mode 

ranked ( l»l,Zo) as calculated with the finite difference 
technique] i.e., <£ t (X,Y) is a finite dimension column vector, 
which contains the values of at all coordinate positions 
(NX, NY) within and on the boundary of the cabin space. 

The normalization of these modes is arbitrary. The maximum value 
achieved at any coordinate position has been chosen to be unity 
and the other values adjusted to retain the computed ratios from 
point-to-point. The normalization is carried into the trans- 
mission prediction with the parameter 


— €*; — 


v ///£> 


In the case where the floor is present the integral is 




cos 






^ * <p * d Art*- 

A MCA 




i-i 
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where j counts over all Interior locations, and 



The program outputs the generalized mass for unit radius, i.e.. 


$6 a)!= E 

J-> 

Iff ‘ 


The volume enclosed is 


V =-a*L 


1T- 


© b +• cos©*, swe, 




giving 


6 




2,[ir— ©o +- Co s © 0 sm 
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3.5 Fuselage Structural Model 

In the previous work [1], the basic structural model was for an 
ideal cylindrical shell, i.e., one without stiffeners or floor. 
In the present work, the inclusion of the orthotropic properties 
of the cylindrical shell and the effects caused by the presence 
of the floor are of concern. 
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Altnougn the level of sophistication achieved in the present 
modeling effort is considerable, it is nevertheless to be con- 
sidered incomplete and further work will be done on the structural 
representation in Phase III. 

The basic analytical technique used to predict the resonance fre- 
quencies of the ring-stringer stiffened shell is that of Mikulas 
and McElnan [7]. In this calculation, the effect of the stiffen- 
ers are considered in terms of a "smeared out" model. This model 
is chosen because it is felt that it gives reasonably good predic- 
tions when applied to a stiffened cylinder of typical aircraft 
construction, i.e., one in which the "density" of rings and 
stringers is sufficient to warrant the "smearing out" or averaging 
of the added stiffness and mass contributions to the shell. The 
Mikulas equation definitely breaks down when the subpanels bet ween 
rings and stringers assume their own modal characteristics. 

The test articles used in the present work are not ideally suited 
to use of the Mikulas formulation because of the sparsity of 
stiffeners. Nevertheless, the formulation is used and considera- 
tion given to its Inaccuracies when applied to the test hardware. 

The second effect of concern is that of the floor. In Ref. [8], it 
is shown that the primary influence of a floor partition in a 
cylindrical shell is to force nodes at the floor line, because of 
the high in-plane stiffness of the floor. Whether the floor is 
attached rigidly or with a pinned connection to the cylinder wall 
is of secondary importance, although some shift in the predicted 
resonance frequencies will occur. Therefore the chosen structural 
model of the fuselage (for Phase II) is a curved, orthotropic 
panel running from floor line-to-f loor line (Figure 8a). The 
calculation is made with the Mikulas equation. 

It is emphasized here that this model is of a preliminary nature. 

A better representation (planned for Phase III) will include the 
restraint offered by the floor itself and the remaining portion of 
the fuselage beneath the floor (see Fig. 8b). Assuming that there 
are two modes which for the configurations of Figs. (8a) and (8b) 
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have the same shape above the floor line, it can be assumed that 
the lower resonance frequency would be that calculated for (8a), 
i.e., without the restraint offered by the floor and fuselage 
below the floor. However, even if the difference in the calcu- 
lated frequencies is disregarded, there are other factors in- 
volved, namely that in the case where the floor and lower fuselage 
participate in the response, the upper fuselage (cabin wall) will 
not respond as well to the acoustic field exciting it because the 
excitation not only has to move the sidewall, in the latter case, 
it also has to move, simultaneously, the floor and fuselage 
beneath the floor. This "restraint" can be built into the model. 
The modes of the more complex configuration are required to do 
this. For this reason, it seems necessary to have a modal analy- 
sis of a configuration at least as complicated as (8b) and this is 
planned for Phase III. One can go another step and add, say, a 
wing. In that case, to shake the cabin wall, the floor, lower 
fuselage and the wing would have to be moved also. A modal analy- 
sis of a model including the wing would be required. In the 
latter case, a finite element modal analysis would be necessary. 

Model for a Curved Orthotropic Panel 

The resonance frequencies are calculated from the equation 


1 4 

u) m L. 

ir 4 D 


M’fH/f + M 4 




H- 


D Xj 



12L 4 (»^) f l+SA s +RA,.+5KAr g 
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where m - total mass/unit area (skin, stringer, and rings smeared out 


Er^r ring bending 

Esls stringer bending 

DA* parameter 

D SLy parameter 

jj> SEfA*' __ ring membrane 

•p: .. _ EjAs Stringer membrane 

parameter 

EVlj£y parameter 

A s = 1 «+ oc 4 {^ s /dT( 1 "+~ p 1 )*' 

A r - H- 8N l (S r /<0 t \~M 

p+i6N 4 (? r /^’-(i+^T 

A= (t+-f’) , --hzp x (|4 V M)^ + 5) 


+ ( i—/**) [ § + s 4 r +- z p- rs ( t +• m)J 


A rs = 4N l o< 1 l-/**')+Z(l+/U)J (Ej/a-) 1, + l6N 4 [l-/* ,t 

+ 2|3‘(i+/*)] (Zy/oC? + 8N l (i-X)(S s /a) 


). 


+ 8N t (l- / M')(? r M + 3ZN 4 0- / H-) 1 (€ r /aXZ s /<0+ l-M 



cabin pressure - exterior pressure 


and H? s are the distance to the centroid from the skin middle 
surface for ring and stringer, respectively , Jl x andJJyare the ring 
and stringer spacings , E,. , E^; X r ,£ s ; A r , A 5 are the elastic moduli!, 1 
moments of inertia, and areas of the stiffeners, and h is the skin thickness 


oC MtT^- / L_ , 

where M is the number of axial halfwaves. 


D— Eln 3 /|^(»-/A. x ) y 


is the skin bending rigidity. 

(3 =2NL/Mira. , 

where N is the number of circumferential half wavelengths. 
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It is again emphasized here that the above equation is strictly 
applicable to a full cylinder. It is also approximately appli- 
cable to an open cylindrical shell having all shear diaphragm 
supports. This boundary condition is obviously not satisfied at 
the cylinder floor Junction. Because of this fact, the theory 
developed by Peterson and Boyd [8] which includes the full cylinder 
with floor, is far superior for present needs and will be inte- 
grated into the model in the future. 

3.6 Calculation of the Structure/Acoustic Coupling Term 
f'(n.r) (Cylinder with Floor ) 

The circumferential structural mode shape for the curved panel 
(floor-to-floor) is given by 

* H (e) = sin = sin , 

where N is the number of circumferential half waves and Ly is the 
length of the panel from floor-to-floor, i.e.. 


L r = (ir-e o '). 


The acoustic mode eigenvector takes on discrete values on the 
boundary (as calculated with the finite difference technique). 
Consider the cylinder boundary point j “(NYjNY). The angle 
is given by 


©j = 


NX-G 

cos0 o -NY-G 
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The acoustic volume mode eigenvector on the boundary line 0j is 

4 >l (Sj 1 ) COS , 

where &< e !») is the ith eigenvector from the finite difference 
calculation in 2 dimensions at boundary point j as. (NX f NV) and 
angle 0j . 


Let 6j^and ©^define the points on the circumference half-way 
between boundary point (NX, hy) and the adjacent boundary points. 
Assume that the eigenvector i) applies over the range 0, • to 


e. 


Then 




Now 




Nit(e-eJ 

z.Ctt-e.’) 


Je 



a. <f> t (6;) f cos .HirC&ij-e.) _ Co NitCe^-e.) Izcn-e,) 

1 L zttr-e.) 2 (ir-e„) J nit 


where is the number of cylinder boundary points. 


- 7 *»- 



Thus , 





Nit(e.j-a) 

icir-e^ 


cos 


zor-e,} J ' 


where, for the case considered (shell length L = cavity length L c » 

r i-cos(Mnnr . i-cos(m-^w 1 
L M+* M-ft J * 


V - 


0 


M-H^ 


If 
noting 


L^U C , ”fa M is given by Eq. 
ing i.ySs ctr-e^ , the abov 

*V> MN)= "rItC ^[ Cos 

J ^ 


(25) of Reference [1]. Finally, 
above reduces to 


HTr(6.i-e») C oC Ntr ^ etiH 
zd r-e e ) 2(t r--‘ ’ 


Lil®s) 1 

*-e 0 ) J 
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where L is the index assigned to the ordered 2D-acoustic mode 
(s is always equal to unity if the floor is present), and M>N are 
the number of axial and circumferential halfwaves of the 
structural mode. 

3.7 Calculation of the Acoustic Lobb Factor 


The calculation of the acoustic loss factor at frequency cO has 
been discussed in Section 3.1 (resulting in Eq.28). This has to 
be formatted properly for the case where the floor is present. 
Also, when the test articles are considered the absorption 
capability of the end caps has to be considered. In both test 
cases where the floor Is present, the end caps are covered with 
fiberglass that Is exposed to the interior. The end cap admit- 
tance calculation is as follows (where the subscript C indicates 
their consideration): 
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In this simple model, the end caps themselves have been considered 
to have a high impedance. The loss factor for the Interior is 


> = 


_ IOC, 


Here 


CO, 


J f ft Jx + 2 


curved 

5ur-f*ct 


end 

Cop 


“ T" ^ [ e *j " e 'j] 

Sof-4«ice J*"’ 


where is the number of boundary points, and 


f ft ^ =^]Z ft JArUi= J2 , 


end 


Are* 


where is the number of Interior points. 
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Therefore at low frequencies 






+ a* $,.(<-)} > (62) 


where the volume mode is rl «= ( 

At high frequencies 

*^\T^ S +z ^ Se J ’ ( 

Sand being the curved surface and end cap areas. 

3.8 Noise Reduction In the Volume Stiffness Controlled Region 


The frequencies of concern here are those confined below the first 
(lowest) resonance frequency of the acoustic modes. Since the 
lengths of the cylinders used in this study exceed their dia- 
meters, the frequencies lie below that where the wavelength in the 
air equals twice the length of the cylinder. For the present 
models, 1.219m (48 inches) and the frequencies lie below about 
140 Hz. In this frequency range, the noise reduction depends on 
the net volume displacement of the shell structure, and on the 
compliance of the cavity (as an airspring). If C e is the compli- 
ance of the cavity (confined air volume) i.e.. 



( 64 ) 
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and 



is the compliance of the shell. 


i. e. , 


c 




( 65 ) 


where Av is the total volume displaced by the shell structures 
[including end-caps] when exterior pressure is applied and 
interior pressure is realized, then the noise reduction is 
found to be 

fog. = 20le 3 [l + ' (66) 

Section 3.8.2 of Ref.[l] may be consulted for the development of 
this result. 

There have been numerous studies [9, 10, 11] of the low frequency 
noise reduction of cylinders (both stiffened and unstiffened) 
where it has been found that the measured noise reductions and 
theoretical predictions do not agree in the extreme low frequency 
range. In fact, in the Phase I study [1] significant over- 
prediction at low frequencies was found to occur. Without 
exception, the (measured) low frequency noise reduction is less 
than the predicted noise reduction. End cap transmission and 
leaks have been postulated as the culprits, but while end caps can 
be significant transmitters, they do not tell the whole story. 
Neither do leaks, which can be prevented. 

The basic question which remains unanswered is why a cylinder's 
noise reduction, typically, turns out to be only 25 or 30 dB when 
theoretically it should be 50 or 60 dB, or more. 
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Now for purposes of the present study, there are two types of 
shells to be considered (i.e., stiffened and unstlff ened) . For 
either, the fundamental modal characteristics lead to the occur- 
rence of the first mode (in flexure) above the lowest acoustic 
mode. That being the case, the compliance of the cylinders is 
dominated by their stiffness. Since the noise reduction is still 
significant, the cylinder compliance is approximately 



( 67 ) 


where is the exterior pressure. Thus to calculate Cp, a 

static analysis might simply be considered. The displacement 
YV,=-W over the cylinder surface is computed due to an applied 
pressure b , and integrated over the surface to obtain Av. 


To begin this study, consider some results for an ideal cylinder. 
The basic partial differential equations for an isotropic cylin- 
drical shell element are given by Kraus [12]. His axial coordinate 
is x (this report has previously used 2: for the axial coordinate). 
The equations of interest are 


d u , I-/4 , \+M . M 
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The shell is assumed to extend from K«0 to X**- 1_ • w>v> and vv 
are the displacements in the X» 9, and ^directions, ^pointing 
outward from the middle surface of the skin, with 

»<«*•(.»•)=•»*(£ +-£t) , 


where 


S zz r 9 . 


is taken to be uniform circumferentially and axially, positive 
Inward. 


The above can be simplified by neglecting the longitudinal and 
circumferential inertia terms, i.e., by setting 

P at- ~ P 'it 1 °‘ 

In the case of a perfectly circular shell , closed on the 

ends, with uniform external pressure, it can be shown [12, p.9^], 
that 

This result also holds true for the case where the radius varies 
slightly. It follows that an approximation can be used 


a* 




and Eq. (70) then reduces to 

Jl. v + w + -i- [ 2 - + 

12 U r ae 


e r 
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Eh 

(71) 
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This differs only slightly from the result one would get if ^u/dx 
had simply been set to zero in Eq.(70); Eq.(7l) is a slightly 
better approximation. The axial strain has little influence on 
the result. 

It also follows that Eq.(69) reduces to 


I-/M , J_ Sv , J_ __ n 


(72) 


Eqs.(71) and (72) are independent of u and thus Eq.(68) is no 
longer of concern. 


The membrane theory solution is obtained by setting V*0 , and 
letting the bending resistance of the shell go to zero (i.e., 
setting the first term in Eq.(71) =0), yielding 

w , i-a«* „ 3V _ i -M/t. h 

Y % E V W Eh ' 


At sufficiently low frequencies, this reduces to the simple form 



trsa. 


Eih 



or 


W — 


EK 



(73) 


The membrane stiffness controlled compliance of the cylinder is 
approximately (neglecting local bending near the ends of the 
cylinder) 
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For the 0.020" thick stiffened cylinder 


Cp = |.2.6» IK 3 / psl 


The compliance of the acoustic cavity is 


c c =— , = 7+i.wy r si 


This gives a noise reduction due to the cylinder wall displacement 
of 

NR ss 20lo 3 j H--^] = SS 4-dB 


As will be seen in Section 4, the actual measurements indicate 
about 30 dB. The noise reduction due to the stiffness compliance 
of the end caps (which will be considered later) varies between 
47 and 49 dB. Adding up the end caps and cylinder membrane 
compliances gives in the range -= 3.9 to 4.6 in 3 /psi or a noise 
reduction of 45.6 - 44.2 dB., which is still well above 30 dB. 


It follows that bending deflections of the cylinder wall have to 
be considered. However, if the radius Y" is constant, W must of 
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necessity be independent of £ and the only bending effects must be 
associated with the axial direction and related to the supporting 
end conditions. The 0 independency of w does not apply when the 
concern is with an orthotropic shell (stiffened with rings and 
stringers). For such a shell, Eq. (71) can be generalized somewhat 
to 



+ 2D, 
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+ CL 
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i av 

tr x 9e . 



0 -M(z) ^ 

' t-x r > 


(75) 


where vn is the average surface mass. Eq.(72) remains unchanged. 


Eq.(75) applies for an orthotropic shell in which the rings and 
stringers have been smeared into the shell. The bending rigidi- 
ties Dy and Og are for the sections perpendicular to the X and & 
(i.e.,s ) axes respectively and Py$is the cross rigidity. 


One solution of Eqs.(72) and (75) is for the case of a free-ended 
cylinder where V is independent of X and 


W 



at every point in the shell. Then 
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In this case, the shell deforms in bending only. The condition 

Sv 

v. = -— , (77) 

is the Rayleigh criterion for inextensible vibrations of a ring. 

If Eq.(77) holds, the middle surface of the shell (I.e., ring) 
does not stretch. One might hypothesize that a shell's end 
conditions could be such as to allow Eq.(77) to hold approximately 
(not perfect supports) but if that is the case, the subpanel 
transmission (i.e., between rings and stringers) would turn out to 
be so great that the predicted noise reduction would be less than 
that measured. 


For instance for a subpanel of the stiffened cylinder used in this 
study, as shown in Figure 9, Eq.(76) would reduce to 

Dv 4 w + m^ ”8) 

which is independent of the curvature effects. The displacement 
in the first mode which satisfies (78) is 
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gives 
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ITS 


[(* p 

~ vn u 



(79) 
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The displaced volume is then 


IT 4 rn aJ'f 


yielding (with the data from Pig. 9) 


( 80 ) 


Cp ^ 77 m*/t si. . 

This leads to a noise reduction 

NR = 2.0 log [ I + - ZO.S'JB 

if only one sub-panel is considered. However, there are 30 sub- 
panels, so the total compliance is 


Cj, = 30(77) = 2310 in3/psi 
giving a noise reduction of only 

NR = 2.0 log [ i + ^~ \ - 2A48. 

If one considers the subpanels clamped along their edges, the 
compliance is 13*64 inVpsi per subpanel leading to an overall 


C.^= 30 x 13.6 - 409 in3/ P si 

and 9 dB. 

Thus throwing out curvature effects leads to an underprediction of 
the noise reduction. Since the inextensibility assumption is 
almost impossible to justify, there are good reasons, as seen 
above, not to invoice it. 
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It follows from all the discussion preceding, that it appears 
membrane stresses are in fact working., but that the deflection of 
the cylinder wall is not obeying the simple result given by 
Eq.(73)» In fact, if the noise reduction is as low as has been 
calculated (less than 10 dB), and the measured is 30 dB, one might 
conclude that membrane stresses have to be present. That being 
the case, it should follow that solutions of the simple equation 


W I-A/2- b 

r* Eh * 


( 81 ) 


should be examined. It also seems crucial to study the effects of 
out-of-roundness (i.e., non-circularity). 


Suppose non-circularity of the cylinder is considered, say by 
letting the radius of the shell be given by 

r - a (i- ef ceV) 3 1 . 


Then at first, the solution appears to be 

w ~ — l ~^ /z 1_ z ^' te) ) ■ 

However, this result is not correct because the parameter Y" in 
Eq.(8l) is the radius of curvature ! 

To distinguish between the radius of curvature and the radius of 
the shell, let the radius of curvature be given by cap R. Then 
Eq.(8l) becomes 


W — l— Al/ 2 - u 

IT'-- ~¥T f 


( 82 ) 
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Also 


L. 

c t =i 0 / 0 -'^ R ‘ J> " ls - 

- L °i r ~ [ RT<xJe • (M) 

En •* q 

[Incidentally, ds^KJe , since 0 is measured from the center of 
the cylinder instead of from the center of curvature (which moves 
about) ] . 

Since the radius of curvature of a shell with a slight variation 
in radius can differ greatly from the radius itself, it can be 
seen from Eq.(83), that if R*^©© , Cp can become very large relative 
to that of the perfect cylinder (thereby reducing the noise 
reduction) . 

In cylindrical coordinates, the radius of curvature is given by 


Now if 



t- - a. (i-6-f(e>). 


(85) 


( 86 ) 


where fee) is a smooth positive function normalized such that 
"f(0) = 1, then 


roa* 
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de* = afePce ' > ’ 

— 

and 

a[ i -efo) +e*(fce) + f 'cl))] 3/z 
{ l +6^f "(6)— Z-f(e)| + <= * [V(e) 4-z f '(e) - 
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■fce)?(e)j j 

or 


O 

N l+-6-f "(e) 

since £'(&> is considered snail. 

(87) 

Suppose now, the out-of-roundness occurs over an arc length of 

shell As=dAO s a0i » as shown In the Figure 10. Then if a length 

* / — 
(axially) is considered, the compliance, C^> for that length is 

according to Eqs.(84) and (87) 

/ r [ &l da 

p - Eh LJ 0 (i+^'Ce))*- 

— 

+/ Jel • 

e, 

(88) 
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This reduces to 


c , = ^p.> [(lir _ ei , +9i] 


= aaL ^ aa -’[»+(6,-e,)]. 


(89) 


where 


€> 


i 



Jo 

(t +-^f "(e)) z 


(90) 


It is obvious from Eq. (89) that whenever the integral ©, exceeds 
the interval & t , that portion (shell length X) transmits better 
than the classical membrane theory predicts. The compliance for 
the whole cylinder becomes. 



Ztt(L-jO(l-7M/Z)0? , f 

Eh ih — L 2ir+ (®- e 0 

2ir(i-v«/z) a. 5 L r. X (©,—©,) l 

Eh L 1 + TT ~W~ J ' (91) 
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If the length of out-of-roundness extends over the entire length 
of the shell, je~L , and 


c p = [ **+(®r ©■)] > 

which is consistent with Eq.(89). 

Before considering evaluation of the integral ©j, it is also 
useful to consider a solution in which the transverse 
displacement v is not assumed to be zero. Prom Eq.(71), the 
membrane result is obtained by setting (alsoOi—O ). 



Eh 


F 


i 


and from Eq.(72), by writing (since R is a function of 0) 


I -vM 
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R deLR\te 




Then substituting Eq.(93) into Eq.(9*0 gives 
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(93) 
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hr a 
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R &©■ L l +■ 

„ 1 — m /2. 

f>. f — afc^'"ce) ] 

EK 

R L (1 +-fef"(e)) a J 


Using Eq.(87), this becomes 

l-'M _ _ (\-/Ufz)€$'\e) 

Z dx z I -t*ef"ce) Eh * 

or 

¥v_ _ _ z-A _P_ e-P'"(e) m 

dx* 1-/4 Eh <4-£-f"te) 


This can be integrated to yield 

vrx,«= + Vj » 

where 

_ _ 2^4 JL e-Ffe) ( (96 

I -A* Eh l+-tf"(S) 

and V 2 , Vi are constants of Integration, (which may be functions 

ore). 

If at some locations (say at the rings) v(x,e)=o , then \/ z and V3 
can be determined. Take 
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Then 


v(o,e) = o 
v u, e) — o . 


V 3 = o 



and one can write 


v(x.,e) = V,(e)-VM 


(97) 


where 


V(x) 

This can be substituted back into Eq.(93) to obtain yy, i«e.. 


w = - 


I-/M/2. . cv2. _ 

Eh r N ae 
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But 


^Vi _ _ 2 .-/M _p_ (H-fef "ce)) fef '(e) - 6.^*(ej 

je i-a» eh ( l ef"(0)) z 


so an approximate solution (neglecting any flexure) is 


Y* 


u-aq e r * z 

Eh(i+. €: f"(e)) a 'L * 



( 98 ) 
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This gives over the length JL 


dd&d'x * 

P 

This integral contains in its denominator the same term 
(n-efW as Eq.(90). If i.e., if 00 over 

any interval, the membrane stress (hoop tension) is broken over 
length Ji . Only the bending resistance remains in the shell in 
this region. This leads to a conclusion that over certain axial 
regions, the cylinder transmission may be totally controlled by 
flexure and over others totally controlled by the hoop tension. 

Now examine the integral for some selected functional forms, 
fee) [refer to Pig. 10]. Suppose 

■f C6) = sin — s’m ^ 0 . 

Then 



IT 

_ a_f dec 

( | - es'm« )*■ 



o o 


where G *■> eir */©;*-. 
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This can be evaluated to yield [13] 


®' = • (99) 
/i ~ 

Since is imaginary for e>l , ©, does not exist in this 

region. This is interpreted as total relief of the membrane 
stress. 


Suppose one chooses 

r . z IT ^ 

+ te)= sin — e • 

This function has zero slope at 0*0 and ©, , which is 
probably a more realistic representation. Then 



where G. — ZfelT/^0^ , giving 
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Fore*. I , the solution is [13] 
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A plot of this function Is shown in Figure 11 


This integral also does not exist (is divergent) for e > \ , thus 

the membrane stress is relieved. 

From the form of these solutions, it is seen that the likelihood 
©| growing large relative to 0 , may be small since combinations 
of 6 (erA) and 0, may not exist which will yield e-cr l . But the 
possibility of £>l may be real (i.e., over the 

range A0=*O +t> & t )• 

Let's choose <& (rather thane) as a criterion for relief of the 
membrane stress. It can be considered an approximate (as opposed 
to an exact) criterion for over an interval A©- . fi is 

plotted versus 0, and 6 in Figure 12. When e> l the membrane 
stress is relieved. 

Now consider some example calculations. Suppose a very slight in- 
dentation of the skin of the 0.020" thick cylinder existed, say 
A ~ O.Ool in. Then according to Figure 12, if 0, 4L 0.0^5 rad 
(2.6°), £>| and the membrane stress is relieved. While 0.001 in. 
may seem small, the amplitude of the membrane displacement W for 
the ideal cylinder is from Eq.(73) 

w = K > 

E.K 

which for the 0.508 m (20 in) diameter cylinder gives 

W *= p * 

At 100 dB, ^ is only 2. 9xl0 -11 psi, giving W =• • In 

comparison A is huge. Unless the pressure is high enough to 
Increase to the same order of magnitude as A , the shell would 
not stretch enough to displace its out-of-roundness . This would 
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apparently require a pressure of about 2 psi or more for the 
.001 in. selected. The interpretation is that very slight 
variations in radius over small angular regions (wrinkles) can be 
highly detrimental at acoustic pressures. 

Now experience tells us that putting dents in a cylinder will not 
change its low frequency noise reduction very much. However in 
our experience, we have never worked with a perfect cylinder. 
Suppose for example that there existed a perfect cylinder with 
wall thickness h — 0.020 in., a radius £*.■= 10 in., and length 
L. 1 ® 48 in. The noise reduction is as has been calculated before, 
i.e. , 55«4 dB. Suppose that over some axial expanse of length X , 
there exists around the circumference a region for which € > 1 . 

For instance, on the stiffened cylinder of this study there is 
about a 2 in. diameter inward depression that is about 0.030 in. 
deep. Setting inches and A = 0.030 in., it is found that 

e,= 4 -= oz - 

4 = — r= O- 003 • 

(X. 

From Fig. 12, e *1.48 > i and thus the membrane stress is 
relieved. The cylinder model is now one where a length L—J?. is 
membrane controlled, the rest (length X) is flexure controlled. 
Consider now a strip running between stiffeners 2 in. wide, and 
from Figure 9 of 11.59 in. length* The compliance for the membrane 
controlled section is 



2ira 3 (l— 4.) 
Eh 


0-/VO 


1.8.1 w\ 3 /psc • 
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The compliance for the strip in flexure between stringers is 


r S _ i 1 !2Afi— /**’) 

=■ 0. 3 1 S' 

Since there are 5 stringers (i.e., 5 strips) 

Cp = s ( 0 . 3 is)= I.S 7 


Also, for the whole cylinder 

Cp= Cp + C f f = 2.78 in 3 /psV , 

and the noise reduction would be 


NR = 20lo<j £ | + = 48.C, dB • 

Thus, this single depression would reduce the noise reduction from 
55.4 dB to 48.6 dB, or about 7 dB, which is a sizable reduction. 


But if the cylinder is not perfect, say without the depression 
its noise reduction is only 35 dB, then its Cp would be (overall) 
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13.45'mYpSc . 


If length 4 was perfectly circular, the portion of the compliance 
above due to it would be 
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With the depression A=0.03 In. and 0 t ®.2, over length , the 
same length's compliance would become 

c.t = i.s7 i^/rsi-- 


Thus with the depression 

13.4* — 0.0S3 = 14-.T7 ’m^/psi. . 


This gives a noise reduction of 3^.1 dB. For the non-ideal 
cylinder, the same depression thus gives only a 0.9 dB reduction 
in the NR. Thus in the non-ideal cylinder case adding a wrinkle 
or depression will not profoundly effect the (already low) noise 
reduction. 


To generalize all of this, a stiffened cylinder with stringers 
will have a compliance given 


C P~ 
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where is length for which e > I • 


Since in general, there may be regions of length with e.^> \ , 
the final expression becomes 
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2lTa 3 ( L. - — /* a) 

Cp= Eh" 
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it,,> fa BaNj-rr 


( 101 ) 
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where =* O if I and inward depressions only are 

considered. 

It should be noted in Eq.(lOl) that larger have more signifi- 

cant impact. Thus a four inch length in which e>| would be far 
more detrimental to the noise reduction than 4 one-inch segments. 


3»9 End Cap Transmission 

The end caps of the test articles are stiffness controlled well 
beyond the fundamental acoustic resonance at 140 Hz. The internal 
pressure due to each end cap's deformation is computed with 



W,(x) dx , 


where £ is an interior point and X lies on the end cap. 

is the interior space Green's function as previously defined. The 

displacement of the end cap under uniform pressure ^ is [14] 


W,(X) =w l M,*=<0=v,(r)= + -r<] , 


where D is the bending rigidity. 


Consider the response -of the first acoustic mode only. Then 
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Substituting, this gives, for both end caps 



. o 4 ft. Co 1 r Tra u l-h+M \1 
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costtz/l^ * 


The space average mean square pressure is obtained with 

< ^>5,*= 77 / f’ctf) f >*(0 Jv » 
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This is averaged over each frequency band of concern, where Jo 
is assumed to be the mean square exterior pressure in the band with 
an approximately uniform frequency content. Then 
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where CO is the band center frequency and Aco— is the 
bandwidth. 


Integration gives 


X. 
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where 


k= ^ , 

"t** * =>J n , 

and 

C0 H = 2TT X 14-0 Ht. * 


>| H is the measured value of the loss factor for the acoustic mode 
(at its resonance (140 Hz)). The noise reduction is obtained 
using 
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3.10 Status of the Interior Noise Program (Phase II) 


A brief status report concerning the computer programs that exist 
for the present (Phase II) calculations is warranted. An overview 
of the software is provided in Figure 13* 

The main program is called FCYLNR and calculates the noise 
reduction of the cylinder with or without floor. If the floor is 
present, it utilizes the output (on tape) from the program CYL2D 
that computes the acoustic modes of the cross section of the 
cylinder with floor. If the floor Is not present, subroutines 
CINDEX and OMEGAN are used as they were for the empty cylinder In 
Phase I. CINDEX Indexes the cavity (cabin) modes. The structural 
modes are computed using the Mikulas and McElman scheme as 
outlined in Section 3.5. There are two programs - MSTRUC which 
prepares input data in proper form to describe the skin, stringer, 
and frame (ring) properties for the stiffened (orthotropic) 
cylinder, and MINCYL which calculates the modes (resonance 
frequencies). Output is again to a tape read by FCYLNR. 

Various subroutines in FCYLNR compute the structure/acoustic 
coupling factors, joint acceptances, radiation loss factors, and 
all other quantities needed In Eqs.(36) and (Ml). These sub- 
routines are structured to .handle the modal data for the cylinder 
with or without floor. Finally, subroutine TRIM computes the trim 
panel transfer matrix data and the required transmission coeffi- 
cients, loss factor, vv^., and the surface admittance Q of the trim. 
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4.0 EXPERIMENTS 


The Phase II study has a two-fold objective: development and 
validation of the analytical models. Section 3 considered the 
development; Section 5 will consider the validation. This section 
considers the experimental program that was undertaken to obtain a 
data base to allow comparisons against predictions for a number of 
different parameters needed in the noise reduction calculation. 
Also the noise reductions of three test articles were measured. 


Most of the experimental work was done at NASA Langley 
Research Center. One simple test was performed at BBN/LA 
to obtain data for checking out the quality of the finite- 
difference calculations of cabin modes when the floor is 
present. The two main tests were: 


1) Measurements of the structural and acoustic loss 
factors of the test cylinders and cavities 

2) Measurement of the noise reduction of the cylinders, 
i.e., the difference between the exterior and 
interior sound pressure levels given a random, 
diffuse (reverberant) exterior acoustic field. 

Other tests done at Langley involved the acquisition of accelera- 
tion response data on the cylinder walls to provide some experi- 
mental insight into the modal characteristics of the shell 
structures. 

4.1 NASA Tests 


The cylinders were tested in a 221m3 (7800 ft3) reverberation 
room having dimensions of about 6x9x4m (20x28x14 ft) with a 
reflection coefficient of about 0.95* The principal room mode was 
near 16 Hz. 
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The testing consisted of the measurement of the noise reductions 
of the test articles (cylinders) and the rate of decay of acoustic 
vibrations in the cylinder cavity and structural vibrations of the 
cylinder wall. 

r 

Description of Models 

The cylinders tested were the unstiffened cylinder with floor and 
the stiffened cylinder with and without floor, with the various 
trims outlined in Section 1. The external dimensions were 0.508 m 
(20 in.) diameter by 1.219 m (48 in) long (see Pigs. 14 and 15). 
The cylinders were equal in weight with one having an unstiffened 
0.0016 m (0.064 in) thick skin of 2024-T3 aluminum. The other had 
the shell weight equally divided between the skin, stringers and 
ring frames. 

The skin was attached to 0.0127 m (0.50 In.) thick annular end 
rings and the 0.41 m (16 in.) diameter opening in the end rings 
was closed with a 0.46 m (18 in.) dia. 0.0127 m (0.5 in.) thick 
cover plate (end cap). The cylinder was airtight except for the 
0.00018 m (0.007 in.) radial clearance between the mike support 
tube and bushings. The bushings and the end caps were attached by 
bolts, the remainder of the attachments and seams were epoxy 
bonds. Both cylinder seams had a 0.00457 m (0.18 in) wide double 
skin thickness; the unstiffened cylinder had a butt-joint with 
outside strap, the stiffened cylinder had a lap joint. 

The stiffened cylinder had 5 stringers, 0.0254 x 0.0254 x 
0.00318 m (1 x 1 x 0.125 in) angles at 72° spacing and 5 ring 
frames, made from 0.00229 m (.09 in) thick sheet, at 0.203 m (8 
In) spacing. The frames were made in two pieces, a 0.034 m (1.34 
in) wide ring with a butt joint, and an annulus with a 0.555 m 
(21.87 in) outside diameter that slipped over the ring to form a 
0.024 x 0.034 x 0.0023 m (0.937 x 1.34 x .09 in) frame cross 
section. 
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In the tests conducted with a floor and trim added to the basic 
cylinder configuration, as shown for the stiffened cylinder in 
Fig. 16, the floor was made from 0.0008 m (0.032 in) thick 
aluminum sheet. The floor was 1.1938 m (47 in) long with flaps on 
the ends to close off the compartment below the floor. A bead of 
silicon rubber (RTV) was used to attach the floor to the skin and 
seal the lower compartment. The floor was stiffened by two 0.0158 
x 0.0158 x 0.00158 m (0.625 x 0.625 x 0.062 in) angles underside 
of the floor. The ends of the stiffeners were bolted to two 
0.0254 x 0.0254 x 0.0032 m (1 x 1 x .125 in) angles on the upper 
side of the floor that provided lateral stiffness and were bolted 
to the end rings to locate the floor about 0.163 m (6.46 in) below 
the cylinder centerline ( 0 O = 49° )• 

The trim insulation was Owens-Corning PF-105 Fiberglas supplied 
0.013 m (0.50 in.) thickness. The supplier. Flight Insulation 
Marietta, Georgia, added a 2 mil vinyl facing on one side. A 
vinyl sound barrier weighting 2.44 kg/m 2 (0.5 lb/ft 2 ) was used. 

The lead-vinyl had a fiberglas cloth facing on one side with a 
total thickness to about 0.00119 m (0.047 in.). 

Instrumentation & Apparatus 

In the noise reduction tests, six microphones were located inside 
the cylinder as shown in Fig. 14. The microphones were mounted on 
a bar attached to a 0.0254 m (1 in) tube on the cylinder center- 
line. The tube could be rotated and translated to position the 
microphone array at the desired measurement location. Exterior 
levels were measured by 2 microphones on floor stands . The 
noise load was supplied by one or two 100 watt speakers (Fig. 17). 


The exterior spectrum was produced by a GenRad noise generator and 
passed through a power amplifier to two Altec "studio monitor" 
speakers. The level was controlled by the volume knob on the 
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noise generator. The amplified microphone output was fed, one 
channel at a time, by a selector switch to an analyzer and digital 
volt meter. Overall levels of all microphones were hand recorded 
from the DVM, one-third octave spectra for selected microphones 
were taped and some NR spectra were plotted. Post test analysis 
and listings of spectra were obtained by recalling the taped 
spectra into the analyzer memory and reading the analyzer memory 
through an IEEE computer interface. 

In the reverberation time measurements, structural and acoustic 
modes were excited by a 5 in diameter speaker mounted on the 
microphone bar. A time history of the signal decay of an interior 
microphone or an accelerometer on the skin was captured by the 
"Hold" function of the analyzer. The time history was then 
recorded by the digital plotter. 

In the mode identification tests, a tone was played through the 
small speaker. One axis of the analyzer scope was connected to 
the speaker and the other to a hand held microphone. The frequen- 
cy was swept until the microphone response peaked. The microphone 
was then moved over the skin surface and the number of slope 
changes in the display per longitudinal and circumferential cir- 
cuit were counted. 

Test Procedure 


The floor and trim were installed in the cylinders. A microphone 
mount short enough to allow the desired rotation was selected, and 
the mikes were installed in the cylinder. The cover plates were 
bolted on the cylinder and the cable access hole sealed with 
modeling clay. The cylinder was suspended from a wood 2 x 4 
"sawhorse" fixture by aircraft bungee cord as shown in Pig. 18. 

The fixture was located off center in the room because the spatial 
gradients in the sound pressure level for the low frequency bands 
were highest near the room center. The exterior microphones were 
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placed about 0.254 m (10 in) from the skin to avoid the increase 
in SPL that occurred for d <£. 0.178m (7 in). Spectra for the 
exterior microphones at the 4 ’'corners'* of the cylinder were 
monitored and minor adjustments to the location and orientation 
■ of’ the cylinder and speakers was made. to obtain approximately 
equal SPL in the 20 and 25 Hz bands. Above these bands the sound 
field was fairly uniform throughout the room. Below these bands 
the interior levels were sometimes lost in the background noise. 

A speaker output level high enough to produce measurable levels on 
the interior microphone nearest the centerline was selected and 
used for all data for the given test configuration. The speaker 
was turned on, OASPL from the DVM was read for all microphones, an 
external spectra was acquired, stored on tape, and transferred to 
analyzer memory "B". An internal spectra was acquired, taped, and 
the difference in current and stored spectra was plotted. After 
comparing the desired microphones the speaker was turned off and 
the microphone array rotated to the next location. 

For reverberation time measurements a 1/3 octave filter, B&K 1614, 
was added between the noise generator and amplifier (Fig. 17) and 
between the microphone selector switch and analyzer. The speaker 
output was down 10 dB for the adjacent 1/3 octave bands on either 
side of the selected one, and down about 50 dB for the second band 
on either side. Time histories of vibration decay were captured 
by simultaneously hitting a switch to disconnect the signal to the 
power amplifier and tripping the ’’Hold” switch on the analyzer. 

If the decay was too irregular to analyze, a tone was substituted 
for the pink noise and the frequency was swept manually until a 
resonance was located. If the decay was smooth the record length 
(time scale factor) and transducer gain was adjusted until the 
decay record filled the CRT screen. The display was checked to 
see that the low-pass filter associated with the time scale knob 
was higher than the frequency of the band being studied, and that 
the decay extended over more time than the 4-cycle time-constant 
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of the filter. The record was then plotted, an arbitrary zero 
time marked on the plot 2 or 3 cycles after the surge due to the 
speaker power switching transient, a curve was hand faired for the 
"average” envelope of the response, and the double amplitude of 
the voltage' scaled off at about 4 selected times, converted to dB 
and plotted on semilog paper. A straight line was faired through 
the data points. Ihe reverberation time, "[^ , i.e. the 60 dB down 
time was calculated from the slope of the faired line. 

Application of trim 

In the stiffened model there was bare metal exposed. The areas 
were the top of the floor, the top of the leg of the stringers, 
and the face of the end rings. The inner side of the cover plate 
had a 0.406 m (16 in.) disk of fiberglass with the facing film 
next to the metal and attached by a few small strips of double- 
back tape. The remainder of the interior surface was covered with 
a double blanket 0.0127 m (0.5 in.) thick made by placing the 
facing of the fiberglass layer against the facing of the lead- 
vinyl layer. The fiberglass side of the blanket was placed next 
to the skin and the lead-vinyl side was placed next to the under- 
side of the floor. The lead-vinyl was cemented to the floor and 
the fiberglass attached by double back tape. For the skin insula- 
tion, the two layers of the blanket were not attached. In the 
lower compartment the blanket was held against the skin by 
gravity. In the upper compartment, small bowed retainer strips, 
covered wagon hoop style, held the blanket against the skin by the 
buckling pressure of the retainer ends against the stringers. The 

retainers compressed the blanket to about 35 "to 70 percent its 
nominal thickness of about 0.014 m (0.55 in.). This left about 
0.01524 m (0.6 in.) of one leg of the stringers bare. Photo- 
graphs showing the exposed stringers are included as Appendix B. 


- 108 - 



4.2 BBN Test 


The sole BBN test consisted of the measurement of the mode shape 
and resonance frequency of an acoustic mode of the cylinder with 
floor. A 0.457 m (18 In.) diameter ceramic cylinder approximately 
0.304 m (12 in) high was used. A 0.019 m (0.75 in.) thick 
fiberboard cover was fitted to the cylinder. Holes were drilled 
for insertion of a 0.0064 m (0.25 in.) microphone. A 0.0064 m 
(0.25 in.) thick plywood partition simulating the floor was taped 
to the side of the cylinder. Excitation of the acoustic mode was 
with a tone, the input being from a small speaker located at the 
top center of the lid of the cylinder. Partitions simulating 
e 0 =0°, 15°, and 30° were tested, for purposes of determining 
the influence of the floor and for comparison to the predictions 
with the finite difference technique. The test considered the 
effect of the floor on the empty cylinder (0,0,2) mode. Results 
and comparisons are presented in the next section. 

4.3 Noise Reduction Data Analysis 

Noise reductions, defined as the difference between the mean 

exterior sound pressure level and the interior level, at measure- 
ment position, i, were calculated and averaged in accordance with 
the relation 

1=1 

where N is the number of sampled (equal volume) interior points. 
The results is the cylinder noise reduction, HR. . The associated 
standard deviation for the noise reduction measurements was also 
calculated with 
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To quantify the accuracy of the estimates, 99# confidence inter- 
vals for the noise reduction estimates were also calculated using 
the relationship 



10 


1 °3io [ 


10 


— NR/10 


VfT ^*"*,0.00$ 


1 


where * s the 0.005 percentage point of the Student "t" 

variable with m = N - 1 degrees-of-f reedom. It was assumed that 
the individual noise reduction estimates were statistically 
Independent, which is believed to be an acceptable assumption. 
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5.0 RESULTS AND COMPARISONS 


In this section, results from various tests are considered in a 
sequence which allows for a gradual confirmation of some of the 
basic analytical models used in the noise reduction predictions 
prior to the comparisons of the predictions with the measured 
noise reductions themselves. Each test article is considered in 
turn. 


5.1 Cavity Modes (Cabin with Floor) 

To begin this study, it is first interesting to examine some basic 
experimental results that were obtained from the test discussed in 
Section 4.2. A very strongly responding acoustic mode in the 
empty ceramic cylinder was identified as the ($,«,s) = (0,0,2) 
mode. The mode shape is [1, p.22] 

<*>(£)= Jo(m 01 r) , 

0 , 0,1 

where the root of (iv^a) = 0 is = 3*83/01, where O. is the 

radius. The theoretical resonance frequency is [l, p. 22] 

XL — ™ - 3.33 C 8 

'0,0,2. OT 3.1T CL 

which for the 0.457 m (18 in) diameter cylinder gives 

fwr * IS “»• 

The measured mode is shown in Figure 19* The resonance frequency 
is almost exactly the theoretical value and the mode shape closely 
matches the theory. 
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In Pig. 20 a small perturbation has been introduced by a 15° floor 
and the same mode has been identified. The resonance frequency 
has shifted upward to 922 Hz and the mode shape shows slight 
bulges appearing at about the 0°, 100° and 260° positions. 

The change in modal pattern is even more pronounced in Pig. 21 
where the results for a 30° floor are shown. The resonance 
frequency has shifted to 930 Hz and the mode shape bears little 
resemblance to that of the no-floor case. 

It was after the results of these experiments were available that 
a decision was made to abandon a perturbation approach to the 
acoustics problem in favor of the finite different technique. 

After the theory of Section 3. 4 was developed, predictions were 
compared against these data to assure the proper working of the 
computer program and the sufficiency of the chosen grid spacing. 

Pig. 22 shows the data of Pig. 21 in a more detailed form. The 
measured resonance frequency is (as previously stated) at 930 Hz. 
However this measurement is for a 0.2286 m (9 in) radius cylinder. 
The test articles are 0.254 m (10 in.) radius. Thus this mode 
would have appeared at (9/1 0) x 930 Hz or 837 Hz if the test had 
been performed in a 0.254 m (10 in) radius cylinder. 

Pig. 23 shows the third symmetric mode for a 0.254 m (10 in) 
radius cylinder calculated with the finite difference technique. 
The computed resonance frequency is 836 Hz. The calculation is 
for a 40° floor. It should have been made for a 30° floor for a 
direct comparison but there was a slight communications problem 
in house. Nevertheless the first twenty symmetric and anti- 
symmetric modes were calculated with the finite difference tech- 
nique, then a search for a symmetric mode closest to the measured 
frequency (837 Hz) was undertaken, and the mode of Pig. 23 was 
found. Comparisons of Pigs. 22 and 23 show that the finite 
difference technique is quite accurate, certainly accurate enough 
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for the present purposes. In fact, the acoustic mode shapes and 
frequencies can be considered as perhaps some of the best data 
that is input to the noise prediction scheme. 

Appendix C contains the first 11 modes computed for a cabin with 
a 0 o = 49° floor per the test articles. The modes are ranked in 
terms of the ”Xj_. Mode 5 (symmetric) is the same mode considered 
in Pigs. 19 through 23- The computed value of the resonance 
frequency for this mode (for Ct = 1 m) is 213 Hz. This gives, for 
a 0.254 m (10 in) radius cylinder, a resonance frequency of 
838 Hz. 

5»2 Noise Reduction Measurements, Predictions, and Comparisons 

As stated previously, the primary purpose of the present -work is 
to demonstrate the feasibility of the calculation of sound trans- 
mission into the interiors of certain complex test articles. The 
noise reduction (in one-third octaves) has been chosen as the 
measure for comparison. The three test articles are: 

1 ) a bare ring-stringer stiffened cylinder without floor or 
trim; wall thickness of 0.000508 m (0.020 inches). 

2) a 0.0016 m (0.063 in.) thick unstiffened cylinder (the 
Phase I test article) modified with a floor partition, lined 
with a 0.0127 m (0-50 in.) thick layer of PF-105 fiberglass 
that is covered with a 0.0000508 m (0.002 in.) vinyl 

film, with a 0.0127 m (0.50 in.) layer of PF-105 fiberglass 
on the inside of each end cap exposed directly to the 
interior. 

3) a 0.000508 m (0.020 in.) ring-stringer stiffened cylinder 
(same as (l) above) with a floor partition lined with a 
simulated trim consisting of a 0.0127 m (0.50 in.) thick 
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layer of PF-105 fiberglass covered with a 0.00119 m 
(0.047 in.) thick layer of lead vinyl weighing 2.44 kg/m2 
(0.50 lb/ft2), with a 0.0127 m (0.50 in.) layer of PF-105 
fiberglass on the inside of each end cap exposed directly to 
the interior. 

In Cases (2) and (3) above the exposed floor surface in the cavity 
(cabin) is bare metal. In Case (3) the (internal) stringers 
protrude somewhat above the trim exposing a height of about 0.1524 
m (0.6 in.) of bare metal on each of five stringers. 

Noise Reduction Measurements 


Fig. 24 shows the equal volume sampling scheme used when the floor 
was present. The sampling was done at 18 positions on three 
measurement planes representing one-fourth of the cylinder volume. 
Because of the nature of the exterior field and the symmetry 
present, the 54 measurement positions are sufficient to compute 
the space average mean square interior level and the ratio of 
exterior level to average interior level, i.e., the noise reduc- 
tion. In reality, this was accomplished by measuring the indivi- 
dual noise reductions nrl at the 54 positions and averaging 
according to the relations given in Section 4*3* The standard 
deviation and the 99# confidence intervals were also based on 
these 54 measurements. In the case of the stiffened cylinder 
without floor, a slightly different measurement procedure was 
followed, with larger sampling volumes, resulting in only 38 
measurement locations. Pigs. 25, 26, and 27 show the measured 
noise reductions. 

5*2.1 Predictions for the Bare 0.020 in. Ring-Stringer 
Stiffened Cylinder 

If insulation (or insulation and lining, i.e., trim) is present, 
the wall conductance can be estimated using Eq.(33) and then the 
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acoustic loss factor can be computed with Eq.(28). Also the added 
structural damping induced by the trim can be determined (ij^). 

The structural loss factor (without trim present) must always be 
measured for input to the computer program. In the case of a bare 
cylinder, the program requires the input of both acoustic and 
structural loss factors. Figs. 28 and 29 show the measured values 
used for the bare ring-stringer stiffened cylinder. 

The acoustic and structural resonance frequencies are also 
required in addition to the loss factors. The acoustic reson- 
ance frequencies for the cylinder without floor are computed with 
the relation 

— C o 9 

where is given on page 22 of Ref.[l]. The structural 

resonance frequencies come from Eq.(60) of this report. Once the 
U) n 's and u) r 's are known, the constants b r » C n * C r > Dnr 
arctan n , arctan^, ln^, and ln r , i.e., everything in the braces 
{ ] in Eq.(36), can be computed. The acoustic mode normalization 
€ n comes from Eqs.(l9) and (20) of Ref.[l]; the modal mass M*. from 
Eq.(23) of Ref.[l]. Since trim is not present ^ | C. W I ^ = 1* 

The acoustic/structure coupling term is given by Eq.(24) 

or (25) of Ref.[l], the joint acceptance J*. , by Eqs. (26), (27) 
and (28) of that same reference. 

Structural Model 


The stiffened cylinder is made up of a thin skin with 5 external 
rings and 5 internal stringers. It is freely supported at its 
ends, and the mode shapes are assumed to be of the form 


, MN. » 

V (e.*) = 


sin 


Mire 

L 


CoS N 0 

sm N0 y 
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Given below is a tabulation of the data used to specify the 
properties of the stiffened cylinder as required by Eq.(60). 


Stiffened Cylinder Properties 
(Dimensions in Meters) 


Cylinder Data 

Length = 1.2192 
Radius = .25375 
Skin Thickness = .508000E-3 
Youngs Modulus = . 72400E+11 N/m2 
Poissons ratio = 0.33 
Density = 2700 kg/m3 
Total mass/unit area = 4.27217 kg/m2 
Stringer Data 

Stringer Spacing = .318867 

Cross-Sectional Area = .151209E-03 

Mass per Unit Length = .408265E+00 

Moment of Inertia About Centroid = .904240E-08 

Distance of Centroid to Skin Middle Surface = -.007768 

Torsional Constant for Stiffener = .46541 6E-09 

Stringer Element Input Data 

Width Height z to Skin Inner surface 

.0254000 .0031750 -.0015875 

.0031750 .0222250 -.0142875 

Frame Data 

Frame Spacing = .203200 

Cross-Sectional Area = .121935E-03 

Mass per Unit Length = .329225E+00 

Moment of Inertia About Centroid = .45371 6E-08 

Distance of Centroid to Skin Middle Surface = -.005207 

Torsional Constant for Stiffener = .200933E-09 

Frame Element Input Data (Frame Elements) 

Width Height z to Skin Inner Surface 

.0342900 .0022860 +.0016510 

.0022860 .0190500 +.0123190 


Resonance frequencies calculated for the "smeared-out" stiffened 
cylinder and an unstiffened 0.02 in. cylinder are shown in 
Figure 30. For the stiffened cylinder the membrane stiffness 
terra dominates until N ^ 5* By the time the bending stiffness 
term dominates, the modes have at least one-half wavelength 
circumferentially between the stringers. As the axial mode number 
increases to MS: 6, the effects of the rings are reduced and the 
subpanel motion begins to dominate. 
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Hu, Gormley, and Lindholm [15] analyzed a ring stiffened cylinder 
using discrete elements for the rings and shell. Their more exact 
analysis suggests that the Mikulas formulation over-predicts the 
resonance frequencies slightly until the sub-panels break-up, at 
which time the Mikulas theory is riot applicable. 

Since the cylinder being analyzed is a thin cylinder (.02 in. 
thick) with relatively heavy stiffeners, the effects of the rings 
dominate for low orders of M and N. The resonance frequencies 
calculated using the ' smeared-out ' rings and stringers are 
therefore used when calculating the cylinder noise reduction, 
since the cavity modes at these frequencies couple to the lower 
order values of M and N. 

Referring to Pig. 30, it is seen that the lowest predicted 
frequency is for the (M,N) = (1,2) mode. Examination of the 
response curve of Pig. 31, made with a shaker attached to the 
skin, shows a resonance frequency at 205 Hz. Another resonance 
occurs at 292 Hz. The first resonance is a structural mode, the 
second is apparently the (2,0,1) acoustic mode (see Table 1, 
p. 84, Ref.[l]). The next significant response is at 407 Hz which 
could be structural or acoustic. On the basis of this measure- 
ment, the (M,N) = (1,2) mode is shifted downward to 205 Hz. The 
(2,2) mode predicted at 660 Hz is shifted to 645 Hz. The (1,3) 
mode is not shifted since it is a more questionable action. 

Noise Reduction Prediction 


Figure 32 gives the calculated noise reduction curve. Also 
plotted are the measured values. Table 1 gives a breakdown of the 
five (5) highest contributing acoustic and structural modal pairs 
in a fashion similar to the manner presented for the unstiffened 
.063 in. cylinder of Phase I. 
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TABLE 1 


Modal Pairs Having Higiest Contributions 
to Interior Level 
Stiffened Cylinder 


Band 
Center 
Freq. (Hz) 

Acoustic 

Mode 

q,n,s Freq. 

Structural 

Mode 

M,N Freq. 

Contribution 

< * , i > 

<p|> 

Overall Noise 

<p?> Reduction 

<p|> nr 

50 

MD CAPS 


.1 122E-04 

.1122E-04 

49-5 

63 

END CAPS 


.1348E-04 

.1348E-04 

48-7 

80 

MD CAPS 


.1949E-04 

.1949E-04 

47.1 

100 

END CAPS 


.4168E-04 

.4168E-04 

43-8 

125 

MD CAPS 


• 5248E-02 

.5248E-02 

22.8 


160 1,0,1 

141.7 

4,0 

3440.6 

.1207E-05 


0,1,1 

396.1 

1,1 

627.4 

•5798E-05 


2,2,1 

715.6 

1,2 

205.0 

.8295E-05 


0,2,1 

657-1 

1,2 

205-0 

•5379E-04 


1,0,1 

141.7 

2,0 

3440.6 

.5640E-04 


END CAPS 



.1219E-01 

.1231E-01 19. 1( 


200 2,1,1 

487.1 

1,1 

627.4 

.2127E-05 


4,2,1 

867.9 

1,2 

205.0 

•9443E-05 


0,1,1 

396.1 

1,1 

627.4 

.2576E-04 


2,2,1 

715.6 

1,2 

205.0 

.5512E-03 


0,2,1 

657.1 

1,2 

205.0 

•3363&-02 


END 

CAPS 



.2512E-05 

•3954E-02 24.03 


- 118 - 






TABLE 1 (continued) 


Modal Pairs Having Highest Contributions 
to Interior Level 
Stiffened Cylinder 


Band 
Center 
Freq. (Hz) 

Acoustic 

Mode 

q,n,s Freq. 

Structural 

Mode* 

M,N Freq. 

Contribution 

<pj> 

<p!> 

Overall 

<p£> 

< P i> 

Noise 

Reduction 

NR 

250 

1,1,1 

420.7 

2,1 

1305.4 

.1571E-05 




2,1,1 

487.1 

1,1 

627.4 

.8492E-05 




2,2,1 

715-6 

1,2 

205.0 

.3074E-04 




0,1,1 

396.1 

1,1 

627.4 

.1271E-03 




0,2,1 

657.1 

1,2 

205.0 

.1902E-03 

.3619E-03 

34.41 


315 

2,0,1 

2,0,1 

2,1,1 

0,2,1 

0,1,1 

283.5 

283-5 

487-1 

657.1 

396.1 

3.0 

1.0 
1,1 
1,2 
1,1 

3440.6 

3440.6 

627.4 

205.0 

627-4 

.2268E-04 

.3606E-04 

.4132E-04 

.9055E-04 

.1265E-02 

. 1 490E-02 

28.27 

400 

0,1,1 

396.1 

3,1 

1784.3 

.202 IE-03 




0,2,1 

657.1 

1,2 

205-0 

.2079E-03 




2,1,1 

487.1 

1,1 

627.4 

.5979E-03 




1,1,1 

420.7 

2,1 

1305.4 

. 1 1 58E-01 




0,1,1 

396.1 

1,1 

627.4 

.6021E+00 

.6148E+00 

2.11 


* N = number of circumferential waves. 
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TABLE 1 (continued) 


Modal Fairs Having Highest Contributions 
to Interior Level 
Stiffened Cylinder 


Band 
Center 
Freq. (Hz) 

Acoustic 

Mode 

q,n,s Freq. 

Structural 

Mode 

M,N Freq. 

Contribution 

<p2 > 

Overall 

<pj> 

Noise 

Reduction 

NR 

<P|> 

<p 2 > 

e 

500 

0,3,1 

903.6 

1,3 

506.5 

•7287E-03 




1,2,1 

672.2 

2,2 

645-0 

•9537E-03 




2,1,1 

487*1 

3,1 

1784-3 

.1981E-02 




0,1,1 

396.1 

1,1 

627-4 

.4825E-02 




2,1,1 

487.1 

1,1 

627.4 

.2987B+00 

•3083E+00 

■ 5.11 


630 

3,2,1 

0,1,1 

4.1.1 
0,2,1 

1.2.1 

782.7 

396.1 
691.6 
657-1 

672.2 

2,2 

1,1 

1,1 

1,2 

2,2 

645-0 

627.4 

627.4 

205.0 

645*0 

.1727E-01 
.2748E-01 
.3063E-O1 
.1787B+00 
. 1 629E+02 

.1658E+02 -12.20 

800 

2,2,1 

715.6 

1,2 

205.0 

.1051E-01 



2,3,1 

947.2 

3,3 

813.9 

.1986E-01 



2,2,1 

715-6 

3,2 

1038.9 

.2099E-01 



4,2,1 

867.9 

3,2 

1038.9 

.3837E-01 



3,2,1 

782.7 

2,2 

645.0 

•9341B-01 

.2140E+00 6.70 
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TABLE 1 (continued) 


Modal Pairs Having Highest Contributions 
to Interior Level 

Stiffened Cylinder 


Band 
Center 
Freq. (Hz) 

Acoustic 

Mode 

q,n,s Freq. 

Structural 

Mode 

M,N Freq. 

Contribution 

<p[> 

< p ‘> 

Overall 

<p|> 

Noise 

Reduction 

NR 

1000 

5,3,1 

1148.5 

4,3 

1085-3 

•3201E-01 




6,2,1 

1074.7 

3,2 

1038.9 

•3448E-01 




0,3,1 

903.8 

3,3 

813. 9 

•3601E-01 




2,3,1 

947.2 

3,3 

813-9 

.1004E+00 




3,3,1 

998.8 

4,3 

1085. 3 

.1738E+00 

•5193&+O0 

2.85 


1250 4,4,1 

1276.8 

5,4 

1370.8 

.6198E-01 


5,3,1 

1148.5 

4,3 

1085-3 

.1198E+00 


3,4,1 

1220.5 

4,4 

1143.2 

.1254&+00 


1,4,1 

1152.7 

4,4 

1143.2 

.2947E+00 


0,4,1 

1380.2 

1,5 

1395-3 

.4214E+00 

. 1 352E+01 -1.31 


At 1600 Hz and above the calculations are made with the high 
frequency result, Eq.(4l). 

Discussion of Results 


Referring to Fig. 32, it is noted that below 125 Hz, in the 
volume stiffness controlled region, a good prediction cannot 
be made because the compliance of the cylinder is determined 
by the relief of the membrane stresses as discussed in Section 
3.8. The predicted values shown in Fig. 32 (below 125 Hz) are 
based on the flexure of the end caps (Section 3*9) • 
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Now suppose there exists a region for which 5 > 1 and the membrane 
stresses are relieved. Then if the length involved is 6. = 0.05 m 
(2 in), the noise reduction at 50 Hz is predicted to be 46 dB, 
i.e., based on the compliance of the shell, Cp , given by Eg. (101). 
If the length is 0.076 m (3 in) the noise reduction is predicted 
to be only 35 dB. For JL = 0.10 m (4 in), the value is only 24 dB. 
However data do not exist that can be used to determine the extent 
of the regions for which § > 1 for the present shell. Thus a 
prediction is really not possible below 125 Hz. 

Above 125 Hz, a fairly good prediction has been made based on the 
NASA measurements of the acoustic and structural loss factors and 
the modal characteristics of the shell. The 630 Hz band is anoma- 
lous due to the predicted presence of a highly resonant (2,2) 
structural mode closely coupled in both wavenumber and frequency 
to a resonant acoustic mode, i.e., the (1,2,1) mode. The (2,2) 
mode has a very high Joint acceptance because it has acoustically 
fast wavespeeds in both the axial and circumferential directions. 
From the data in Table 1, the contribution of the (2,2) and 
(1,2,1) modes can be seen to significantly dominate the predic- 
tion. The noise reduction can easily be computed without their 
participation. The result is the dashed curve in Figure 32 which 
falls much closer to the measurement. If the prediction is viewed 
with the (2,2) mode contribution suppressed, it is seen that the 
prediction follows the measurement with reasonable accuracy from 
125 Hz out. The noise reduction climbs to a maximum in the 250 Hz 
band where there are no resonant acoustic nor structural modes, 
then nose-dives to a value of only 2 dB at 400 Hz. This is 
followed by a plateau where the noise reduction remains below 10 
dB before increasing slightly in the last few bands. 

5-2.2 Predictions for the 0.063 in. Unstiffened Cylinder with 
Floor and Insulation 

As in the previous case, the acoustic and structural loss factors 
are required; also the resonance frequencies, acoustic and 
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structural mode normalizations, coupling factors , and joint 

acceptances. Since the floor is present, the finite difference 
acoustic modal data are used and since insulation is included, the 
trim transfer coefficient is needed. 

Figure 33 shows the calculated structural loss factor for two 
different values of the polyester film loss factor that appears 
in the calculation of through the parameter C w , given by the 
trim transfer matrix coefficients on p. 34* Also shown is the 
measured value of that was used in Phase I. The Phase I data 
were available and were used as the bare cylinder input. However, 
since the insulation was present, and the cylinder was so lightly 
damped without it, the value of V could have been set to zero 
without affecting the calculated value of significantly. 
Generally there is an overprediction in >|p , however, the fact 
that the predicted and measured values generally lie in the range 
between 0.01 and 0.10 is encouraging, that is, basically the right 
order of magnitude is being predicted. 


The acoustic loss factor measurements and predictions are shown in 
Fig*34« As can be seen the amount of damping afforded by the 
flexure of the 0.00005 m (0.002 in) vinyl film is quite drama- 
tic. The calculation for various values of y |t are shown. The 
equations used to calculate were those of Section 3*7, specifi- 
cally Eqs.(62) and (63) where the conductance came from Eq.(33)« 

A good comparison was never achieved and as a result, the envelope 
of the measured acoustic loss factor indicated by the horizontal 
line at = 0.05 and sloping off above 2000 Hz was used as compu- 
ter input. As can be seen, the measured acoustic loss factors 
bounce around and identifying the few data points with any speci- 
fic acoustic modes was impossible. 
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Figure 35 shows the computed value of the trim transfer 
coefficient (i.e., -10 log^) . Again the influence of the value 
of is apparent. For the present work, a value of = 0.5 was 
selected to prevent a large negative excursion in the range 
between 1000 and 2000 Hz. 

The acoustic modes were computed with Eg. (58) and the structural 
modes with Eq.(60), where the model was a curved isotropic panel 
running from floor line to floor line. The acoustic mode normal- 
ization was determined with Eq.(59) and the coupling factors with 
Eq.(6l). The joint acceptances came from Eq.(26) of Ref.[l] 
where j^is computed with Eq.(26) as is, and with Eq.(26) 
changed such that N replaces M and Ly replaces Lx • As before 

.1 ‘2 .1 

J f J MN Jm’Jm* 

Figure 36 shows the predicted noise reduction and Table 2 presents 
the computations in detail. 
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TABLE 2 


Modal Bairs Having Highest Contributions 
to Interior Level 

Unstiffened Cylinder with Floor and Insulation 


Band 

Center 

Freq. 

Acoustic 

Mode 

q,i (or n)Freq. 

Structural 

Mode* 

M,N Freq. 

Contribution 

< P 2> 

<p2> 

*e 

Overall Noise 

<p?> Reduction 

<F> ® 

50 

END CAPS 


.1072E-04 

.1072E-04 

49-7 

65 

END CAPS 


.1288E-04 

.1288E-04 

48.9 

80 

END CAPS 


.1905E-O4 

.1905E-04 

47.2 

100 

END CAPS 


•4074E-04 

•4074E-04 

45-9 

125 

END CAPS 


.2630E-01 

. 2650E-01 

15.8 


160 0,5 

857.5 

1,5 

140.9 

.5008E-04 


0,5 

657-5 

1,4 

187.5 

.6510E-04 


2,0 

285.5 

1,5 

140.9 

.71 25E-04 


0,2 

444.9 

1,5 

140.9 

.8867E-04 


0,1 

579.6 

1,4 

187-7 

.9752E-04 


END CAPS 



. 1 049E-02 

.1761E-02 27-5 


200 2,1 

475-8 

1,4 

187.7 

.9104E-04 


2,5 

716.0 

1,4 

187.7 

.1166E-05 


0,6 

957.1 

1,4 

187.7 

.2262E-05 


0,5 

657-5 

1,4 

167.7 

•7076E-05 


0,1 

579.6 

1,4 

187.7 

.1 165E-02 


END 

CAPS 



.251 IE-04 

.2555E-02 25.9 


* N = number of circumferential half-waves. 
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TABLE 2 (continued) 


Nodal Pairs Having Highest Contributions 
to Interior Level 

Unstiffened Cylinder with Floor and Insulation 


Band 

Center 

Freq. 

Acoustic 

Mode 

q,i (or n)Freq. 

Structural 

Mode 

M,N Freq. 

Contribution 

< P 2> 

Overall 

<p2> 

Noise 

Reduction 

NR 



250 

0,1 

379.6 

1,2 

620.1 

.5523 E -04 




0,1 

379-6 

1,4 

187.7 

.8284 E -04 




2,0 

285.5 

1,9 

245-1 

.11 33 B -03 




0,2 

444-9 

1,3 

309.6 

.1549 B -03 




2,0 

285.5 

1,3 

309.6 

.1709 E -02 

.2466 E -02 

26.08 


315 

2,2 

527-5 

1,3 

309.6 

.8007 E -03 



0,1 

379-6 

1,2 

620.1 

.8769 E -03 



0,4 

703.7 

1,3 

309.6 

.2269 E -02 



0,2 

444.9 

1,3 

309-6 

.1094 E -01 



2,0 

283-5 

1,3 

309.6 

•3384 E -01 

•5033 B -01 12.98 


400 0,1 

379.6 

1,4 

187-7 

.6424 E -03 


3,0 

425.2 

2,5 

434.4 

.301 2 E -02 


1,2 

466.9 

2,5 

434.4 

•3064 E -02 


0,2 

444-9 

1,3 

309.6 

.1165 E -01 


0,1 

379.6 

1,2 

620.1 

.1 576 E -01 

•3719 E -01 14.30 
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TABLE 2 (continued) 


Nodal Fairs Having Highest Contributions 
to Interior Level 

Unstiffened Cylinder with Floor and Insulation 


Band 

Center 

Freq. 

Acoustic 

Mode 

q,i (or n)Freq. 

Structural 

Mode 

M,N Freq. 

Contribution 

<p|> 

<pi> 

Overall 

< P 2> 

<p|> 

Noise 

Reduction 

NR 

500 

2,2 

527-5 

1,3 

309-6 

.1337E-02 




0,1 

579.6 

1,2 

620.1 

.200 IE-02 




1,2 

466.9 

2,5 

434-4 

•3485B-02 




2,1 

473-8 

1,2 

620.1 

.695 IE-02 




0,2 

444.9 

1,3 

309.6 

•9709E-02 

.2751 E-01 

15-60 


630 

3,1 

570.0 

2,4 

622.6 

•3651E-02 



2,1 

473-8 

1,2 

620.1 

.3664E-02 



0,1 

379-6 

1,2 

620.1 

.6530E-02 



1,3 

672.6 

2,4 

622.6 

.5600E-01 



0,3 

657.5 

1,2 

620.1 

•7078E-01 

.1552F+00 8.09 


800 3,3 

783-0 

2,4 

622.6 

.1863E-02 


1,4 

717-9 

2,3 

960.2 

.2184E-02 


3,4 

822.2 

2,3 

960.2 

.2477E-02 


0,5 

837.5 

3,5 

829.4 

.4922E-02 


2,5 

884.1 

3,5 

829.4 

.1144E-01 

•3643E-01 14.39 
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TABLE 2 (continued) 


Modal Fairs Having Highest Contributions 
to Interior Level 

Unstiffened Cylinder with Floor and Insulation 


r 

Band 

Center 

Preq. 

Acoustic 

Mode 

q,i (or n)Freq. 

Structural 

Mode 

M,N Freq. 

Contribution 

< P 2> 

Overall 

<p|> 

< Pe > 

Noise 

Reduction 

NR 

1000 

2,7 

1001.3 

3,5 

829.4 

.2361 E-02 




3,6 

1029-0 

4,6 

971.3 

.2678E-02 


- 


2,6 

979.0 

3,4 

1128.2 

•3087E-02 




2,5 

884.1 

3,5 

829.4 

•359 IE-02 




4,6 

1095-2 

3,4 

1128.2 

.6439E-02 

•4339E-01 

13-83 


1250 

5,9 

1322.7 

6,7 

1351.8 

.2767E-02 



3,9 

1195.0 

4,5 

1228.3 

•3792E-02 



5,7 

1193-5 

4,5 

1228.3 

. 5530E-02 



4,6 

1095-2 

3,4 

1128.2 

•5S38E-02 



4,10 

1329.8 

5,6 

1297.4 

.6185E-02 

.5306B-01 12.75 


Discussion of Results 


In general the results of the predictions are good, the shapes 
of the measured and predicted curves being quite similar. Below 
125 Hz the predictions are based on the end caps. The plateau 
appearing in the measurement between 250 and 1000 Hz has been 
predicted. This is followed by a rise in the noise reduction and 
then a sharp increase between 3150 and 5000 Hz due to the 2 mil 

vinyl film covering the fiberglass on the cylinder wall. 

This conclusion is evident from examination of Fig. 35. In the 
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160 and 200 Hz bands, the measured noise reduction exceeds the 
predicted value, however, the sharp drop at 250 Hz is basically 
predicted once the density of resonant acoustic modes is 
sufficient. Another interesting observation is that the 
prediction generally tries to follow the measurement beyond 1000 
Hz, that is rising slightly to a value above 20 dB before the 
effect of the vinyl film is felt. It is emphasized that this 
calculation is made using the complex modal characteristics of the 
cabin space created by the floor partition. To our knowledge no 
sound transmission prediction of comparable complexity has been 
attempted before. 

5.2.3 Predictions for the 0.02 in Ring-Stringer Stiffened 
Cylinder with Floor and Trim 

The calculation of the acoustic modes was made with Eq.(58). The 
finite difference data were used because of the presence of the 
floor. The structural modes were computed with Eq.(60), where the 
model was that of a curved orthotropic panel running from floor 
line to floor line. The acoustic mode normalization was deter- 
mined with Eq.(59) and the coupling factors with Eq. .(61). The 
joint acceptances came from Eq.(26) of Ref.[l] and were calcu- 
lated in the same manner as they were previously for the 0.065 in 
cylinder with floor (see 5*2.2). 

Acoustic Loss Factors and Trim Transfer Coefficient 


Predictions were made of the acoustic loss factors for various 
values of nt of the lead-vinyl covering the fiberglass on the 
sidewall. These are compared to measurements in Figure 57- If 
there is assumed to be no dissipation in the trim ( = 0), the 
predicted acoustic loss factors do not match well with the 
measured data. However by assuming the trim panel is dissipative 
and taking a loss factor of 0.5 or greater, fairly good predic- 
tions can be achieved. Not only are the wide excursions in the 
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predictions of the acoustic loss factor suppressed, but also those 
present in the trim transfer coefficient as seen in Figure 38. A 
value of *u = 0.5 was chosen for the noise reduction computation. 

Structural Loss Factors 


Figure 39 shows the measured structural loss factors. Also shown 
are the predicted values of Since the trim was present when 

the structural loss factor was measured, theoretically r\£, was 
measured. However, the excitation was with a speaker inside the 
cylinder and ideally, to determine the exterior should be 
excited allowing the trim to react to that response to provide the 
added damping. Although the predicted «V is very high at 400 Hz 
and below, it is seen from Fig. 38 that this is precisely the 
region over which -10 log / f^ is negative, that is, when an 
increase in the power flow to the interior due to the presence of 
the trim is predicted (as compared to the no trim case) . The 
terms dominating the are determined by the value of the para- 
meter Cy, which becomes large and acts to reduce the power flow, 
i.e., in opposition to the effects caused by being greater than 
unity. At frequencies where is less than unity (i.e., 

- 10 log >0), tj£_ drops by at least an order of magnitude. 

Influence of Stringer Exposure 

At high frequencies, when the cavity response is resonant and the 
structure response is also, the exposed stringers (protruding 
above the trim) can be a significant contributor to the interior 
sound level. In fact, if the power flow through the trim itself 
becomes sufficiently low, the stringers can become the primary 
contributing source. In the present case, the power flow through 
the trim is predicted to be very low at high frequencies leading 
to high noise reductions. Since the primary contributing struc- 
tural modes are all acoustically fast, the radiation off the 
stringers is simply the ratio of exposed area of the stringers to 
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the transmitting area (with trim). Thus the acoustic power inflow 
due to radiation from the stringers is estimated to be (approxi- 
mately) given by the power flow computed to pass through the 
cylinder wall had the trim not been present, multiplied by the 
ratio of the exposed stringer area to the cylinder wall area 
(floor-to-floor) , but allowing for the increased structural damp- 
ing present (i.e., Vr ) being created by the trim. This leads to 



■trim, on 

AsWiMTS x W,„ W " . 

Aeyl.wall I 


The dissipation in the cavity is almost exactly the same in both 
cases, i.e., with or without stringer exposure. Thus the space 
average mean square pressure in the interior due to stringer 
exposure is given by 


StyWiacrs 

<p^r a _ 

no exposure 

r 5,* 


sfcruuitrs 

w-.u 


> 


or the difference in interior sound pressure levels is 


|0 | 0 . A^rs 
Acyl, vwll 


ioi^icH 2 . 


Thus the interior level with stringers exposed is estimated to be 
greater than the interior level without exposure by an amount in 
dB equal to 




trs 

y 


Figure 40 shows the predicted noise reductions and Tables 3 and 4 
present the computations in detail. 
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TABLE 3 


Nodal Fairs Having Highest Contributions 
to Interior Level 

Stiffened Cylinder with Floor and Trim 


Band 

Center 

Freq. 

Acoustic 

Mode 

q,i (or n)Freq. 

Structural 

Mode* 

M,N Freq. 

Contribution 

<P|> 

<p 2 > 

*e 

Overall 

<p|> 

<p z > 

*e 

Noise 

Reduction 

NR 

50 

HID CAPS 


.1445E-04 

.1445E-04 

48.4 

63 

END CAPS 


.1737E-04 

.1737E-04 

47.6 

80 

END CAPS 


.2512E-04 

.2512E-04 

46.0 

100 

HID CAPS 


•537QE-04 

.5370E-04 

42.7 

125 

END CAPS 


.169SE-02 

. 1 698E-02 

27.7 


160 2,0 

283.5 

1,1 

961.7 

.1854E-05 


1,0 

141-7 

2,5 

676.7 

•3578E-05 


0,1 

380.0 

1,2 

416.5 

.1236E-04 


1,0 

141.7 

2,3 

636.6 

.3208E-04 


1,0 

141.7 

2,1 

1727-7 

•7993E-04 


HID 

CAPS 



.1698E-02 

.1828E-02 27-3 


200 0,2 

445.3 

1,3 

298.8 

.41 78E-05 


1,0 

141.7 

2,1 

1727.7 

•5409E-O5 


2,1 

474.1 

1,2 

416.5 

.6116E-05 


2,0 

283.5 

1,1 

961.7 

.1117B-04 


0,1 

380.0 

1,2 

416.5 

.81 78E-04 

.1276E-03 38.9' 


* N = number of circumferential half-waves. 
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TABLE 3 (continued) 


Modal Pairs Having High est Contributions 
to Interior Level 

Stiffened Cylinder with Floor and Trim 


Band 

Center 

Freq. 

Acoustic 

Mode 

q,i (or n)Freq. 

Structural 

Mode 

M,N Freq. 

Contribution 
<P l> 

w 

Overall 

<p|> 

W 

Noise 

Reduction 

NR 

250 

0,2 

445-3 

1,3 

298.8 

.2687 E -04 




2,1 

474.1 

1,2 

416.5 

.411 IE -04 




2,0 

283-5 

1,3 

298.8 

•3309 E -03 




2,0 

283-5 

1,1 

961.7 

.7044 E -03 




0,1 

380.0 

1,2 

416.5 

•7146 E -03 

.1914 E -02 

27.18 


315 

0,2 

445-3 

1,3 

298.8 

.2222 E -03 



2,1 

474.1 

1,2 

416.5 

.2415 E -03 



2,0 

283.5 

1,1 

961.7 

. 1 357 E -02 



2,0 

283-5 

1,3 

298.8 

. 1 400 E -02 



0,1 

380.0 

1,2 

416.5 

.1 161 E -01 

. 1 536 E -01 18.14 



467-3 

2,3 

636.6 

.1862 E -03 


1,1 

405-5 

2,2 

967-0 

.2452 E -03 


2,1 

474.1 

1,2 

416.5 

.7846 E -03 


0,2 

445-3 

1,3 

298.8 

.1253 E -02 


0,1 

380.0 

1,2 

416.5 

.7657 E -02 

.1081 E -01 19-66 
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TABLE 3 (continued) 


Modal Pairs Having Highest Contributions 
to Interior Level 

Stiffened Cylinder with Floor and Trim 


Band 

Center 

Freq. 

Acoustic 

Mode 

q,i (or n)Freq. 

Structural 

Mode 

M,N Freq. 

Contribution 

cp|> 

<PI> 

Overall Noise 

<p|> Reduction 

<p|> ® 

500 

2,2 

527.9 

1,3 

298.0 

.2160E-03 



0,1 

380.0 

1,2 

416.5 

•4274E-03 



0,2 

445-3 

1,3 

298.8 

•7272E-03 



1.2 

467-3 

2,3 

636.6 

.8114B-03 



2,1 

474.1 

1,2 

416.5 

.1321E-02 

•4200E-02 23.77 


630 

0,3 

0,3 

1.3 

1.4 
3,2 

658.1 

658.1 

673-2 

718.5 

615-7 

1,2 

1.4 

2.4 
2,3 
2,3 

416.5 

397.4 

552.4 

636.6 
636.6 

.301 IE-03 
•3685E-03 
.6084E-03 
.8336E-03 
.1105E-G2 

•4402E-02 

23-56 

800 

0,5 

838.3 

3,5 

808.5 

.1508E-O3 




2,3 

716.8 

3,4 

81 1.9 

.2682E-03 




2,5 

884-9 

3,5 

808.5 

•3294E-03 




4,3 

868.6 

3,4 

811.9 

.4294E-03 




1,4 

71 8. 5 

2,3 

636.6 

.6685E-03 

.2874E-02 

25.42 
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TABLE 3 (continued) 


Modal Fairs Having Highest Contributions 
to Interior Level 

Stiffened Cylinder with Floor and Trim 


Band 

Center 

Freq. 

Acoustic 

Mode 

q,i (or n)Freq. 

Structural 

Mode 

M,N Freq. 

Contribution 

<p|> 

<5? 

Overall 

< P 2> 

W 

Noise 

Reduction 

NR 

1000 

4,4 

904.2 

3,3 

1004.2 

.51733-04 




2,6 

979-9 

3,6 

960.1 

.6108E-04 




5,5 

1097.7 

4,5 

1032.0 

.62533-04 




3,6 

1029.9 

4,4 

1108.5 

.8596E-04 




3,7 

1051.1 

4,5 

1032.0 

.2210E-03 

. 1 250E-02 

29.03 


1250 

4,9 

1253.4 

3,7 

1218.7 

. 1 426E-04 



3,9 

1196.0 

2,7 

1174-3 

.1446E-04 



1,9 

1126.8 

2,7 

1174.3 

.2042E-04 



6,7 

1263.4 

5,5 

1309.0 

.27243-04 



0,9 

1117.9 

1,7 

1165.1 

.38163-04 

.34673-03 34.60 
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TABLE 4 

Effect of Stringer Exposure* 


, X 

Frequency 

(Hz) 

Computed NR 
w/o Stringer 
Exposure 

-10 log|C w | 2 

Estimate of NR with 
Stringer Exposure + 

1000 

29*0 

20.8 

17.0 

1250 

34*6 

25.4 

18.0 

1600 

40.8 

29.6 

20.0 

2000 

50.9 

33.8 

25-9 

2500 

59.9 

37-5 

31-2 

3150 

65-3 

40.7 

33-4 

4000 

71.6 

43.5 

36.9 

5000 

74.7 

45-9 

37.6 


p 

+ Stringers exposed area = 0.1*15 m 

2 

Cylinder wall area (trim covered) - 1.09 m 
= 0,5 

Discussion of Results 


As in the previous cases, the prediction below 125 Hz is based on 
the end caps' transmission. While the prediction below 125 Hz is 
not particularly of concern, it can be seen that the measured 
noise reduction is in closer agreement to the prediction than it 
was for the bare stiffened cylinder. It can be seen from Table 5 
that at 125 and 160 Hz the end caps dominate. This is due to 
their intimate coupling with the (1,0) acoustic mode at 141*7 Hz. 
Incidentally, recall the (1,0) acoustic mode is q = 1 , i = 0 
where i is the assigned order of the 2 dimensional modal pattern 
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in the cylinder cross section when the floor is present (mode 0 
of Appendix C). By 200 Hz, the end cap contribution can be 
ignored. The measured and predicted noise reductions rise 
sharply due to the absense of any resonant acoustic modes and 
thgn fall dramatically as has been observed in the previous two 
comparisons. However the curves rapidly bend upward because the 
trim lining is so heavy (see Pig. 38). The prediction is that 
the noise reduction could rise to over 70 dB by 5000 Hz if the 
stringers were not exposed. Actually, it is a pretty remote 
possibility that it would ever achieve that level even if the 
stringers were not exposed because of other shorting paths 
(vibration of the end caps, exposed metal such as the floor, 
etc.) In the present case, the stringers are obviously the 
shorting mechanism (Table 4)* The transition between the two 
curves in Figure 40 (stringers exposed versus not exposed) is 
not clear, although it could probably have been determined; it 
did not seem important to do so. It is felt that the present 
interpretation of the results is reasonably accurate, and that 
basically all the structure in the curves has been explained 
with the present model. 

5*3 Statiatical Analysis of Prediction Error 

The remaining question to be answered is whether there is a bias 
in the predictions, that is, does the analytical model (with the 
best input data that can be generated) tend to over or under 
predict the noise reduction on the average. 

Table 5 contains the measured and predicted noise reductions for 
the three test articles. It includes the noise reductions for 
frequencies at and above 125 Hz, and excludes the transition 
region data for the stiffened cylinder with floor and trim where 
the stringers' exposure first begins to contribute significantly. 
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Of interest is whether there is a statistically significant 
difference, on the average, between the predictions and the 
measurements for these data, i.e., are the predictions biased? 


TABLE 5 

Predicted vermis Measured Noise Reductions 


Band .02 Stiffened .06 Unstiffened .02 Stiffened 

Center w/o floor w/floor w/ floor & Trim 

Freq. Predicted Measured Predicted Measured Predicted Measured 


125 

22.8 

22.7 

29.8 

27-9 

27.7 

30.4 

160 

19.1 

19.2 

27.5 

34.9 

27.7 

31.9 

200 

24.0 

18.3 

25-9 

35.1 

38.9 

36.7 

250 

34-4 

24-9 

26.1 

17.6 

27.2 

30.9 

315 

28.3 

11.5 

13-0 

16.6 

18.1 

27.0 

400 

2.1 

2.6 

14.3 

15-4 

19-7 

20.0 

500 

5.1 

5-3 

15-6 

16.4 

23-8 

16.9 

650 

-12.2 

4.2 

8.1 

17-8 

23.6 

21.3 

800 

6.7 

6.3 

14.4 

16.2 

25.4 

22.1 

1000 

2.8 

6.1 

13*6 

15-0 

Transition 

region 

1250 

~1-3 

3*9 

12.8 

21.5 

data are not 

1600 

4.2 

4.1 

19-3 

23-3 

included . 


2000 

5-7 

4.6 

20.7 

23-5 

25.9 

30.4 

2500 

7.3 

3.8 

24.2 

24-1 

31.2 

32.1 

3150 

8.2 

3-3 

24.0 

22.1 

33-4 

31.6 

4000 

10.6 

8.0 

31.1 

29.2 

36.9 

36.7 

5000 

12.3 

10.0 

34.4 

36.3 

37.6 

39.4 
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To determine if there is any bias, the differences ; i = 1 » 2 , 
...» n, between the predictions and measurements in each 
one-third octave band are computed and their mean and standard 
deviation determined by pooling all data from the three tests, 
providing a sample size of n = 48. The sample mean and standard 
deviation are 

A =-0.56<te • s -=5.46 <16. 

Under the hypothesis that there is no discrepancy on the 
average, that is, 




O, 


the sampling distribution of A becomes 



> 


where is the student "t" variable with n-1 = 47 degrees of 

freedom. The acceptance region for the hypothesis above is 
given by 


YrT < R 

For a two-sided test at the o< = 0.05 level of significance, 

"k 47 * 0015 ' = 2.01: soS must fall between -1.59 and +1.59 dB. 
Since the sample mean 2i (=- 0.56 dB) falls within the 
acceptance region, there is no reason to question the validity 
of the hypothesis, i.e., there is not a statistically 
significant difference between the predictions and 'measurements 
on the average (at the cC = 0.05 level of significance). 
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On the other hand, there is a substantial random error in the 
predictions versus measurements as indicated by the standard 
deviation of s = 5*5 dB for the A* t values. This means that one 
can expect a discrepancy of more than 5 or 6 dB for about 
one-third of the predictions. 

In the future, output from the interior noise program should 
include, in addition to the predicted noise level, a measure of 
the quality of the prediction based upon a statistical analysis 
of the comparison data accumulated in the various tests. This 
measure would be the confidence intervals, which might be 
computed for the various frequency bands by scaling-up the test 
articles to the fuselage diameter of concern. It is felt that 
there will continue to be a substantial random error although 
its existence will almost always be associated with input data 
deficiencies. 
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APPENDIX A 


Transfer Matrix 
for Trim Insulation 
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Transfer Matrix for Trim Insulation 


Consider a harmonic pressure wave moving from left to right in 
the insulation of Figure 1 (imagine that the skin and lining are 
not present). When the wave reaches x = L, it is reflected. If 
the pressure amplitude in the incident wave at x = L is ^ and 
the pressure amplitude in the reflected wave is then the 
pressure at any point x measured from left to right is (dropping 
the dependence) 


f>Cx)= 


i(k+l«)(x-L) , 

+ fV 


_L(k+\.«0( x-L) 

e , 


(A1) 


where k is the acoustic wavenumber in the insulation and <X is 
the decay constant. Note, that when x = L, 


Let k . Then 




e 


+ 



~V(L.— x) 

e 


(A2) 


The particle velocity is 


Y(X) - 



V(u- x) 
€. 



-r(u-x) 

c > 


(A3) 


where W is the wave impedance in the material, that is 


-1HH- 



w = p/v , 


where p is the pressure and v is the particle velocity for a 
wave running in the material when the material has unlimited 
extent (W is independent of x) . 

Then, when x = L, 


p 

vW = - 

'c- Fv . 

\N 

Now let 


f>(0) = f, 

V(o) = v, 


II 

3 

> 

Also define the impedance 

as 



Then, 


But 



1 + 


\ 

2 Z +W / 


(A4) 



h 



+ 


Zi+w 
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Therefore 


2a.fi 

2* cosh yL +\NsmhYL 



(A5) 


Also, it is found that 


Let 



Wv, 

SmhVL +* W/ cosh Yl- 


(A6) 


Cosh VL_ — C 
si^vh L— 3 . 


Using Eq^s. (A4), (A5), and (A6), after some algebra, it is found 
that 



ws 

c 2 -s z 


V, 


t 


and 



S/w 


h 


c*-s* 



so 




C - \NS 

to 

/ 

-S/W c 



This is the desired transfer matrix in terms of input and output 
pressures and velocities. 
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APPENDIX B 


Photographs of Test Articles 
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Stiffened Cylinder in Test Rig . 
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APPENDIX C 


Finite Difference Results 
Acoustic Modal Patterns and Resonance Frequencies 
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MODE 6. ANTISYMMETRIC 
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MODE I 


SYMMETRIC 
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Exterior 


FIGURE 1. 



TRIM MODEL 


INSULATION AND LINING 


- 166 - 




Glass Fiber Batts 
Total Thickness: 10.2 cm 

(9.61 kg/m^) 


FIGURE 2. TYPICAL WIDE BODY TRIM INSTALLATION 
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FIGURE 3. 



Frequency, Hi 


ACOUSTICAL PARAMETERS FOR OWENS-CORN I NG 
FIBERGLAS TYPE PF 105. FIBER DIAMETER 
IS 1 MICRON AND BULK DENSITY IS 9.6 kg/m 3 
(0.6 lb/ft 3 ) [4] 
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100 1000 10,000 
Frequency (Hz) 


FIGURE 4. COMPARISON OF TRIM TRANSFER COEFFICIENT 
PREDICTION WITH MEASUREMENT 
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FIGURE 6 




FINITE DIFFERENCE GRID NOMENCLATURE 
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FIGURE 7. GRID COORDINATE SPECIFICATION 
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(b) 


FIGURE 8. FUSELAGE STRUCTURAL MODELS 
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FIGURE 9 


SUBPANEL 

jf x = 0.169m (6.6 in) 
ly = 0.294m (11-6 in) 
h = 0.0005m (0.02 in) 


CYLINDER SUBPANEL dimensions 


(SKIN ONLY) 
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FIGURE 10. CYLINDER IMPERFECTION MODEL 
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Membrane Stress 
is Rel ieved Totally 
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FIGURE 11. RATIO OF THE INTEGRAL 0. TO INTERVAL 6 
VERSUS e 1 
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FIGURE 13. STATUS OF CYLINDER NOISE REDUCTION PROGRAM (PHASE II) 
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FIGURE 14. UNSTIFFENED CYLINDER - MODEL DETAILS AND MICROPHONE LOCATIONS 
(DIMENSIONS IN CENTIMETERS) 
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STIFFENED CYLINDER WITH FLOOR 


GR 1382 


Altec 9846-8A 
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FIGURE 17. NOISE REDUCTION TEST INSTRUMENTATION 
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180° 



FIGURE 19. MODAL PATTERN FOR THE (0,0,2) MODE 
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180° 



FIGURE 21. MODAL PATTERN FOR A 30° FLOOR 
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SKETCH OF 3 CYLINDER MEASUREMENT PLANES 
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SECTION A-A. MEASUREMENT LOCATIONS AT EACH 
OF 3 AXIAL SECTIONS 


FIGURE 24. EQUAL VOLUME SAMPLING SCHEME 
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NOISE REDUCTION, DD 


60 


STIFFENED .02 IN CYLINDER WITHOUT FLOOR OR TRIM 
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FIGURE 25. MEASURED NOISE REDUCTION OF THE BARE 0.020 IN. 
RING-STRINGER STIFFENED CYLINDER 


- 190 - 









NOISE REDUCTION, DB 


UNSTIFFENED .06 IN CYLINDER WITH FLOOR AND TRIM 



FIGURE 26. MEASURED NOISE REDUCTION OF THE 0.063 IN. 
UNSTIFFENED CYLINDER WITH FLOOR AND 
INSULATION 
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STIFFENED .02 IN CYLINDER UITH FLOOR AND TRIM 



FIGURE 27. MEASURED NOISE REDUCTION OF THE 0.020 IN. 

RING-STRINGER STIFFENED CYLINDER WITH 
FLOOR AND TRIM 
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Structural Loss Factor 


MEASUREMENTS 

▲ Mike Near Skin (Tone) 

• Accelerometer on Skin (1/3 Octave) 

Computer Input Based on Measured 
Values of Tr (Mike Near Skin Only) 
Averaged Over the Frequency Band 
(4 = Estimated) 



Frequency, Hz 


STRUCTURAL LOSS FACTORS OF THE STIFFENED CYLINDER 
WITHOUT FLOOR 
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Acoustic Loss Factor 



20 100 1000 10,000 
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FIGURE 29. ACOUSTIC LOSS FACTORS OF THE STIFFENED CYLINDER 
WITHOUT FLOOR 
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Frequency, Hz 












ONE-THIRD OCTAVE BAND CENTER FREQUENCY, HZ 


FIGURE 32. COMPARISON OF PREDICTED AND MEASURED NOISE 
REDUCTIONS, 0.02 IN. STIFFENED CYLINDER 
WITHOUT FLOOR OR TRIM 
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Structural Loss Factor, 
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FIGURE 33. STRUCTURAL LOSS FACTORS OF THE UNSTIFFENED CYLINDER 
WITH FLOOR AND INSULATION 
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Acoustic Loss Factor 



FIGURE 34. ACOUSTIC LOSS FACTORS OF THE UNSTIFFENED CYLINDER 
WITH FLOOR AND INSULATION 
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UNSTIFFENED .06 IN CYL. UITH FLOOR AND INSULATION 
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FIGURE 36. COMPARISON OF PREDICTED AND MEASURED NOISE 
REDUCTIONS, 0.063 IN. UNSTIFFENED CYLINDER 
WITH FLOOR AND INSULATION 
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FIGURE 38. TRANSFER COEFFICIENT FOR THE STIFFENED CYLINDER 
WITH FLOOR AND TRIM 
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Structural Loss Factor 



Frequency, Hz 


FIGURE 39. STRUCTURAL LOSS FACTORS FOR THE STIFFENED CYLINDER 
WITH FLOOR AND TRIM 
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NOISE REDUCTION, DB 



ONE-THIRD OCTAVE BAND CENTER FREQUENCY, HZ 


FIGURE 40. COMPARISON OF PREDICTED AND MEASURED NOISE 
REDUCTIONS, 0.020 IN. STIFFENED CYLINDER 
WITH FLOOR AND TRIM (STRINGERS EXPOSED) 
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